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NEW SOLUTIONS OF THE EQUATION OF 
RADIATIVE TRANSFER 


LYMAN SPITZER, JR. 


ABSTRACT 

Solutions of the equation of radiative transfer are given for two assumed variations 
of , with optical depth, valid for n, increasing or decreasing with 7. The first appears 
as a generalization of the Eddington solution. The second avoids the negative absorp- 
tion inherent in previous investigations and leads to a simple formula convenient for 
numerical applications. For certain values of the constants the formulae are expressible 
in terms of simple powers and exponents. 

Discussion of the equation of radiative transfer as applied to the 
formation of a spectral line has hitherto been limited almost entirely 
to the case in which 7,, the ratio of the line and the continuous ab- 
sorption coefficients, does not change with optical depth. Approxi- 
mate treatments for small variations have recently been given by 
B. Stroémgren,' and extensive numerical integrations for variable 7, 
have been carried out by Pannekoek? and others. The only rigorous 


solution, however, has been given by Eddington in the special case 


where 7 is, as usual, the optical depth in the adjacent continuous 
spectrum. 

It is of interest to note that rigorous solutions are available under 
a considerable variety of assumptions as to the functional form of 


Ap. J., 86, 1, 1937. 
2M.N., 91, 139, 1930; ibid., p. 519, 1931. 3 Ibid., 89, 620, 1929. 
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n. The formulae assumed for yn, contain in each case three param- 
eters, and thus the results are applicable to a wide range of 
physical situations. With these solutions it is now possible to make 
a direct mathematical investigation of the effect on absorption lines 
of any simple monotonic variation of the absorption coefficients. We 
give these solutions here for reference as they are apparently not 
available elsewhere in the literature. 

The equation of radiative transfer may be written, following 
Strémgren, 


_ 3(1 + em) | (1) 


dt? I+ I+ 


where ¢,, the optical depth in the line, is defined by 
dt, = (1 + nv)dt = -—(1 + nv) (2) 


We assume 


1 + 
BUt) = av + polo, 


and make the substitution 
y(t) J (ty) (a, + 


If we assume also that 7, is a function of ¢,, the resulting differ- 
ential equation for y may be solved in two general cases. We shall 
neglect ey, in (1), as when this is important (1 + en,)/(1 + 7.) be- 
comes simply ¢, usually attributed to ionization cycles and hence 
roughly constant with ¢,. We then set, for the first solution (solution 
I), 

1+ m = A(t + Dr)’. (3) 


Substituting (3) in (2) and integrating for ¢,, (1) finally becomes 


oy —s/(s+1) 
(1 + mt,) y, (4) 
where 


m=(s+1) 4. 
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Equation (4) may be solved by the substitutions 


as 2 [3 1/2y 
w 


provided s is not equal to —1 or —2, and we find 
y = + oK,(w)} , 


where J,(w) and A,(w) are the usual Bessel functions of imaginary 
argument. Consideration of the boundary conditions J, = 2H, at 
the surface and y equal o at infinity determines c, and c, and leads 
finally to the formula 


m 
where 
2V3A (s+ 1)D 
for D(s + 2) >o, 
= 
VF 
— 7) for D(s + 2) <o. 


In the simplified case that B, is given* by a,(1 + $7), 2p,/3a, is 
simply d7/dt,, evaluated at an appropriate point. Taking this point 
to be t,=1, an “effective depth” of the line, we find 


dt, A A 
4 This is the case contemplated in the definition (5) of r,. For other situations we 
must replace a, by the intensity in the continuous spectrum just outside the line, and 
formulae (5), (8), (9), (12), and (13) became multiplied by the factor 
V3 
2 


I+ 


2 


V3 ( 
34, 
where 8, is the continuous spectrum analogue of /,. 


s+ 2’ 
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A consideration of the analysis as s approaches —1 and « shows 
that (5) and (6) are valid in each limit. In the latter case (s + 1)D 
must approach some constant K and for (3) and (6) we have 


The former case corresponds to the one given by Eddington. The 
difference between his formula and (5) arises from his assumption 
that B, is linear in 7 instead of, as here, in /,. 

For s equal to —4, —4,-—$, —%, —'f,...., etc., 7, is express- 
ible in terms of elementary functions, as y becomes half-integral. 
For D(s + 2) > oas well, F,(¢) becomes a rational function of ¢; in 


the simplest case, s = — 4, we find for positive D, using (7), 
6V34 + 2D I (8) 
oA oA , 


the approximation holding provided D is less than A and A is large. 
Formula (8) corresponds closely to the physical situation represented 
by Eddington’s solution (s = —1) and is considerably simpler to 
apply. For the limiting case s = — 2, (5) no longer holds, but (4) 
becomes a homogeneous linear equation which can be solved direct- 


ly, giving 


(9) 


where the plus sign goes with positive D and vice versa; 2p,/3a, 


may be expressed by (7) as before. 
To determine A, D, and s in practical situations, we may use the 


values of n, and its derivatives with respect to 7, taken at some ap- 


| 

| 
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propriate point. Taking this point to be r = o, and dropping the 
subscript v, we find the formulae 


A=1+ n(o) 
[n'(o) 
[n’(o) — + nfo)]’ (10) 
_ (o) 


where and n’’(r) stand for dy,/dr and d?n,/dr?, respectively. 
The second general solution (II) corresponds to the assumption 


=L+ (11) 


Making the substitution 
w= V3M (u/2) ty 


we find 
y = als(w) + .Ka(w) , 


where 5 = 2V 3L/ u, and J;(w) and K;(w) are again Bessel functions 
of imaginary argument and of the first and second kinds, respective- 
ly. With boundary conditions as before, we find 


ape\t | 3 ) 
(0. 


= = 
M 
rt cs 


where, for u positive, 


(12) 


T;(w) 


G3(w) = 


For M negative both V M and G; become imaginary, and in (12) we 
can take the absolute value of M and for G; the corresponding ratio 
of J;(w) and its derivative. In either case, when 6 is equal to a 


| 
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half-integer, G;(w) becomes expressible in terms of powers, and 
either trigonometric or hyperbolic functions of w. 
Thus, for 6 equal to 3, we have for positive M 


I + 
= (13) 
1+ 3 


2( via coth v1) 


The 2p,/3a, factor, expressed in terms of L, M, and u by the method 
of (7), becomes simply L + Me“, and for 7,(o) large the entire 
second term in the numerator may be neglected. For negative M 
we have merely to replace ‘“‘coth”’ by “cot.” 

There is in this case no simple accurate expression for 7, as an 
explicit function of r. This is not important, however, since we pro- 
ceed as in (10), determining L, M, and u from the derivatives of 7, 
with respect to 7, evaluated at r = o. As before we find 


3ln’(o)P — + n(0)] 
+ n(o)}? 


n’ (0) 
at (14) 
I 


In applying formulae (5) and (12) the physical meaning of as- 
sumptions (3) and (11), respectively, must be kept clearly in mind. 
A must, of course, be greater than 1 if 7,(0) is to exceed zero; the 
values of D and s are not, however, so restricted. With assumption 
I for 1 + », we distinguish four particular cases. 

1. D>o,s>o0. mv increases with 7, and 7 and /, go to infinity together. There 
are no singularities, and we may use (5) without hesitation. The condition on s 
gives, from (10), 9’’(0) < [n’(o)}?/[1 + n(0)]. 

2. D<o, s<o. Here m increases with r but becomes infinite for 7 = —1/D; 
beyond this point it becomes negative, or imaginary, etc., depending upon s. 
But if s is less than —1 (case 2a), ty becomes infinite at this point and thus 7 


| 

| | 

| 

M = | 
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beyond becomes irrelevant and (5) is applicable. For s > — 1 (case 26) this 
is not true, and, in fact, the boundary conditions cannot even be satisfied. In 


2a we have n’"(o) > 2 [n’(o)}?/[1 + n(0)]. 

3. D>o, s<o. mv decreases with r. It approaches —1 eventually, but for D 
sufficiently small this is of no importance. For D large compared to A, however, 
the temperature gradient at t,=1 becomes very large, and we have an emission 
line; i.e., the negative absorption has been brought too near the surface. In 
general, this case must be used with suspicion, unless one is willing to fit solu- 
tions at t = (A-*/* — 1)/D, as in Eddington’s original paper. Note that for 
s<—1, ty remains finite as t approaches infinity. 

4. D<o,s>o. As in case 2 for s> —1, this case has nonphysical peculiarities 
and cannot be used, except in the limiting exponential case s= © , for which the 
singularity disappears and the boundary conditions may be satisfied. This suf- 
fers, however, from the same objection as case 3, since 7» approaches —1 with 
increasing T. 

Thus, we see that (5) may be used with assurance only in cases 1 
and 2a when 7, is an increasing function of the optical depth. Simi- 
lar cases may be drawn up for solution II and formula (12). Cor- 
responding to the restriction A >1 above, we have here the general 
condition 1 > L+M > o. It is clear that when uw is negative, n, 
approaches — 1, in one case decreasing monotonically and in the other 
becoming infinite for finite 7. The formula for Gs(w) would be differ- 
ent here, but since this type of variation has already bee.: discussed 
in connection with (5), and since it is of limited physical interest in 
any event, we shall not pursue it farther. 

For positive u we note that if 7,() is to be positive, L must lie 
between o and 1. Through (14) this together with the condition that 
u be positive imposes restrictions on (0). For n’(o) positive we 
find < 2[n’(o)}?/[1 + n(o)]; for n’(o) negative, > [3 + 
1/n(o)][n’(o) ?/[1 + n(o)]. When M is negative the condition of pos- 
itive M + L implies that the absolute value of M must be less than 
that of L. Thus the argument of Js(5VM /L) is less than its order, 
and J;(w) is a monotonic function in the range considered. 

We see that for n’(o) negative we have a valid solution only for 
n’’(o) greater than the critical value given above. For n’(o) positive, 
however, we have not only solution II but also I, cases 1 and 2a, as 
well, and can thus find ¢, for any values of n’’(0). For (0) <[n’(o) ?/ 
[1 + n(o)] both formula (5) and formula (12) are equally valid. In 
general, both will presumably lead to very nearly the same result for 


| 
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given n(0), n’(o), and 7’’(o), but this conclusion is difficult to check. 
For y and 6 are given by different combinations of these quantities, 
and the only comparison possible is a purely numerical rather than 
an analytical one. In application, moreover, approximations must be 
made, since the Bessel functions are tabulated for only certain orders; 
either the values of the constants must be slightly changed to yield 
tabulated functions or interpolation must be resorted to. At the 
worst a good approximation may be had using (5) and (12) as two 
parameter formulae, taking n(0) and 7’(o) as given, and adjusting 

n’’(o) to give the most easily calculated results. Applying this pro- 
pea to (13), for instance, we set 6=} or u= =4V3L in (14), and 
find a simple cubic equation for VL in terms of n(o) and n’(o). For 

n’(o) negative this has a single positive solution, but for y’(o) positive 
there may be two such or even none. M is found directly from VL, 
n(o),and n’(o). The same method is easily extended to any value of 6. 


PRINCETON, N.J. 
September 10, 1937 
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UNIDENTIFIED INTERSTELLAR LINES IN THE 
YELLOW AND RED* 


PAUL W. MERRILL AND O. C. WILSON 


ABSTRACT 


Six unidentified interstellar lines (Table 1) are discussed. Following a brief history 
of previous observations, Tables 2 and 3 give widths and central absorptions of three 
of the lines. The total intensity of \ 6284, which, in the mean, equals the average for 
the detached lines D1 and D2, shows a somewhat stronger correlation with stellar 
color excess than does the intensity of the D lines. The widths of the unidentified lines, 
which are rather diffuse and not sharp like the D lines, and other facts make an atomic 
origin improbable. The lines may be portions of molecular bands, as yet unidentified. 
Several facts, however, suggest that the lines are produced by small solid particles, per- 
haps closely related to those that cause space reddening. 


INTRODUCTION 


In 1920, Miss Mary L. Heger at Lick Observatory measured two 
faint absorption lines, \ 5780 and X 5797, in several B-type spectra, 
notably that of ¢ Persei.t She evidently recognized the possibility 
that they might be “‘stationary”’ but considered the data inconclu- 
sive. About the same time W. H. Wright recorded a weak line at 
d 6283 in the spectrum of a Cygni? but marked it “probably of at- 
mospheric origin.” 

More recently the lines AX 5780, 5797, 6284, and 6614 were ob- 
served at Mount Wilson, in the spectrum of HD 183143,' but, on 
account of their lack of sharpness, their interstellar origin was not 
immediately recognized. Their occurrence in numerous spectra of 
types Oa~A4 later roused suspicion, however, and correlations of 
their intensities with distance’ and color excess were then soon 
noted. 

The nonstellar origin of \ 5780 and A 6284 was definitely shown by 
their failure to participate in the periodic displacements of lines of 

* Contributions from the Mount Wilson Observatory, Carnegie Institution of Wash- 
ington, No. 582. 

t Lick Obs. Bull., 10, 146, 1921. 

2 Ibid., p. 109. 

3 Mt. W. Contr., No. 409; Ap. J., 72, 98, 1930. 4 Pub. A.S.P., 46, 206, 1934. 
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the spectroscopic binary Boss 61425 and by the fact that the observed 
wave lengths appear, in general, to be independent of the motions of 
the stars in whose spectra the lines occur.° 


WAVE LENGTHS 


Although the positions of these lines have been measured on many 
spectrograms, the accuracy of the mean wave lengths is still not 
very high. On half the plates the lines are near the limit of visibility, 
and, even with spectrograms of good quality where the lines are fair- 
ly intense, the probable error of a single plate (two independent 
measures) is about o.2 A, the slightly diffuse character of the lines 
preventing greater precision. Except for the high contrast and fine 
grain of the III C emulsion,’ with which most of the negatives were 
made, the accuracy would have been still lower. The wave lengths 
computed on the assumption that the displacements are the same as 
those of the D lines in the same spectra are in Table 1. 


TABLE 1 


WAVE LENGTHS OF INTERSTELLAR LINES 


r Intensity No. Plates 
0202500 I— 9 
620600: 104... I 14 
6983-092 6 219 


The lines \ 6203 and \ 6270, which were not included in the orig- 
inal investigation, occur on spectrograms on which d 6284 is strong. 
Reduction of the measures on the assumption that the lines are of 
stellar origin gives the wave lengths 6202.98 + 0.12 and 6270.10 + 
0.08. Since the probable errors obtained by using velocities derived 
from the interstellar D lines (see Table 1) are only half these values, 
the evidence for interstellar origin is strong. 


5 Mt. W. Contr., No. 536; Ap. J., 83, 126, 1936. 
6 Pub. A.S.P., 48, 179, 1936. 7 Supplied by the Eastman Kodak Co. 
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LINE STRUCTURE 

The interstellar D lines, like H and K, are extremely sharp and 
well defined, their apparent shapes undoubtedly being determined 
largely by the limitations of the spectrograph and of the photo- 
graphic plate. The unidentified lines under discussion are, however, 
quite different in character, being wider and softer, with diffuse 
edges. They resemble stellar lines of intermediate sharpness, which, 


45780 A5797 


D3 D2 Di 


| 


OzgAtm. A6284 


Fic. 1.—Tracings of interstellar lines in the spectrum of HD 183143 


although not extremely broad and hazy, might be described as 
“somewhat diffuse” or ‘‘not sharply defined”’ (see Fig. 1; see also 
Pl. X and Fig. 2 of Mt. W. Contr., No. 576).8 

The actual widths of the lines, as well as the total absorptions or 
equivalent widths, were determined from photometric tracings 
(scale, 100 times original, or 0.336 A/mm) of a large number of pho- 
tographs by noting the points at which the line profiles begin to drop 


Ap. J., 86, 274, 1937. 
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below the level of the adjacent continuous spectrum. Although line 
edges cannot be located with precision on individual spectrograms, 
the observations are sufficiently numerous to yield good mean val- 


TABLE 2 
MEASURED WIDTH IN ANGSTROMS 
EW 5780 5707 6284 D2,1 
4.6 (3.8) 6.9 2.2 
A 
12.0 
8.0 - = 
< 
Dp 
0.0 
0.0 0.4 0.8 1.2A 


Equivalent width 


Fic. 2.—Actual width and equivalent width. The dispersion of the original spec- 
trograms is 33.6 A/mm. 


ues. To secure homogeneity all the photometric measurements were 
made by the same observer (O. C. W.). 

The results are shown in Table 2 and in Figure 2. For a given 
equivalent width, \ 6284 is wider than \ 5780, while \ 5797 is prob- 
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ably narrower. The measured widths of the D lines, included for 
comparison, depend to a large degree upon the purity of the spec- 
trum and the properties of the photographic plate, as well as on the 
intrinsic line widths. 


TABLE 3 
CENTRAL ABSORPTION IN PERCENTAGE 

EW 5780 5707 6284 D2,1 Siu 
5 6 4 9 5 


Central absorption 


0.0 0.4 0.8 ¥.2 1.6A 
Equivalent width 


Fic. 3.—Central absorption and equivalent width 


A similar comparison of the measured central absorptions appears 
in Table 3 and in Figure 3. The values for the D lines are undoubted- 
ly too small because of limitations imposed by the spectrograph and 
the photographic emulsion. The values for the other lines, however, 
should be reasonable approximations to the intrinsic central ab- 


| | 
| 
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sorptions.? Since the number of images measured is considerable 
and the scatter of individual values is not large, the means, and 
especially their relative values, should be reliable. For a given equiv- 
alent width, 6284 is shallower than \ 5780, and X 5797 is deeper, 
while all three lines are markedly shallower than the D lines. These 
results are in harmony with those found from the actual widths."° 
The stellar Sz 1 lines, \ 6347 and \ 6371, included in Table 3 for 
comparison, are decidedly deeper than \ 5780 and \ 6284. The 
structure of \ 5797, on the other hand, as far as indicated by the 
meager data, does not differ radically from that of the silicon lines. 
Little can be said concerning the structure of \ 6203 and \ 6270; 
they are perhaps slightly more diffuse than \ 5797. 

The lack of identical structure in AA 5780, 5797, and 6284 indi- 
cates that, if these lines belong to an atomic multiplet, it is one with 
quite unusual properties. 

Concerning the symmetry of the lines no very definite statement 
can be made. \ 5780 usually appears symmetrical, although in about 
ten instances an apparent asymmetry is noticeable on the tracings, 
and in all but one of these the violet edge is steeper, suggesting that 
possibly a faint companion lies toward the red from the main com- 
ponent. We are inclined to think, however, that as a first approxi- 
mation the line is to be regarded as symmetrical. \ 5797 is too weak 
and narrow to give a positive result and \ 6284 is too involved with 
near-by lines of the alpha band of terrestrial oxygen. 


TOTAL INTENSITY 


The integrated intensities of the lines \X 5780, 5797, and 6284, 
expressed in equivalent widths or angstroms of complete absorption, 
are recorded for individual stars in Mt. W. Contr., No. 576.5 A 
graphical comparison in Figure 4 shows that throughout the ob- 
served range the intensity of \ 5780 is, on the average, very nearly 
one-half that of \ 6284 and slightly more than three times that of 


9 Unless, of course, the lines are formed by several narrow, unresolved components. 


‘0 The agreement is not a mathematical necessity because a line with a narrow core 
might conceivably be both wider and deeper than another one of the same equivalent 
width. 


" Largely from spectra of c stars. 
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5797. The intensity of \ 6614, for which only one or two measure- 
ments are available, appears to be between that of \ 5780 and of 
5797. The scatter of the relative intensities, although consider- 
able, does not certainly exceed errors of measurement, which, be- 
cause of lack of sharpness in the lines, are considerable. The inten- 
sity of \ 5797 is below that required for the best determinations, 


> 


° 
> 


Equivalent width \ 5780 


0.0 0.4 1.2 1.6A 
Equivalent width A 5797 and A 6284 


Fic. 4.—Comparison of equivalent widths of \ 5780, \ 5797, and A 6284 


TABLE 4 
TOTAL INTENSITIES OF \ 5797 AND A 6270 


EA 
STAR HD 

5797 » 6270 
4841 0.16 
.09 
20041 .14 
183143 .18 +33 
BOSS GY 223385 0.12 0.10 


while measurements of \ 6284 are made difficult by the presence, in 
the violet wing, of lines belonging to the head of the alpha band of 
terrestrial oxygen. 

According to eye estimates, the intensity of \ 6270 is nearly equal 
to that of \ 5797, while that of \ 6203 is somewhat less. Total in- 
tensities measured from photometric tracings of a few spectra are 


in Table 4. 
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In the spectrum of HD 183143 (see Fig. 1) the intensity of \ 5797 
appears unusually low relative to that of the other unidentified de- 
tached lines. 

Since the intensities of \ 5780 and \ 6284 are closely correlated, a 
combination value has been used in the following discussions. A 
smooth curve, drawn through the plot (Fig. 4) representing ob- 


A 


1.6 


~ 
N 


Equivalent width \ 6284 
> 


0.0 0.4 0.8 1.2A 
Equivalent width 0.5 (D2+D1) 
Fic. 5.—Comparison of equivalent widths of \ 6284 and the D lines 


served pairs of values, was used to obtain from each measurement of 
5780 a hypothetical intensity of. \ 6284 which was combined 
(weight one-half) with the measured intensity. The result is de- 
noted by J’ \ 6284. 

Upon plotting the equivalent width J’) 6284 against that of the 
D lines, 0.5 (D2 + Dr), called J’D (Fig. 5), a strong correlation ap- 
pears, although the scatter is obviously greater than that caused by 
error of measurement. The relatively great intensity of \ 6284, 
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equal in the mean to that of the average of the two D lines, may 
cause surprise. It would perhaps not be anticipated from visual in- 
spection of the negatives, for the eye often underestimates the total 
intensity of a broad, shallow line compared with a narrow, deep one. 
The correlation coefficient for /’D and J’) 6284, from 127 stars, is 
+0.71 + 0.03. 

A statistical dependence of the intensity of \ 6284 upon stellar 
distance (length of light-path), to be expected from the facts now 


| 

N 
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| 
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Distance 


Fic. 6.—Stellar distance and intensity of \ 6284 


available, is shown by Figure 6. The observed irregularities are not 
markedly greater than those exhibited by the calcium” or sodium’ 
lines, although the measurements are probably less accurate. For 
all the lines the scatter at the greater distances and among the more 
intense lines is large. Because of the relatively small central absorp- 
tion, the equivalent widths of \ 6284 in various spectra should be 
nearly in direct proportion to the total numbers of absorbing par- 
ticles encountered by the light-beam."! The scatter in the distance- 

2 Mt. W. Contr., No. 573; Ap. J., 86, 136, 1937. 

13 Mt. W. Contr., No. 569; Ap. J., 86, 28, 1937. 

4 Provided the line is not formed by narrow, unresolved components. 
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intensity relationship would thus measure differences in the average 
concentration of particles in different directions. With the narrow 


calcium or sodium lines, on the other hand, variations in velocity 


along the light-path play an important role, and the relationship 
becomes more complicated. 


LINE INTENSITY AND COLOR EXCESS 


Because interstellar line intensity and space reddening both in- 
crease with distance, they should exhibit a positive correlation with 


T 
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Fic. 7.—Color excess and intensity of the D lines 


each other. This expectation has been verified for the K line by 
E. G. Williams's and by R. F. Sanford."® 

Eighty-five stars with observed intensities of \ 6284 and the D 
lines have colors measured by Stebbins and Huffer’ or by W. Beck- 
er.‘ Scatter diagrams for color excess against J’D and against 
I’) 6284 are shown in Figures 7 and 8, respectively, in both of which 
a positive correlation is obvious. 

1s Mt. W. Contr., No. 487; Ap. J., 79, 280, 1934. 

© Mt. W. Contr., No. 573; Ap. J., 86, 136, 1937. 

17 Pub. Washburn Obs., 15, 217, 1934. Messrs. Stebbins and Huffer have kindly sent 
us additional unpublished data. 

8 Veriffentlichungen der Universitdtssternwarte zu Berlin-Babelsberg, 10, 1, 1933. 
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Because the identification of \ 6284 is as yet uncertain, it is im- 
portant to follow all clues which may bear on the origin of the line. 
Thus it is desirable to determine whether the intensity of \ 6284 or 
that of the D lines is more closely correlated with color excess (space 
reddening). The appearance of Figures 7 and 8 favors a higher cor- 
relation for \ 6284—a conclusion borne out by the correlation co- 


0.4 
3 
*| rs ee 
ee 
| 
0.0 0.4 0.8 1.2 1.6A 


Equivalent width \ 6284 


Fic. 8.—Color excess and intensity of \ 6284 


efficients given in Table 5, which have been computed from the 
standard formula, 


For the 85 stars with complete data the coefficient for color excess 
and I’ \ 6284 is +0.81 + 0.03, while that for color excess and I’D 
is only +0.61 + 0.05. Similar results are found when the 85 stars 
are divided into two groups with /’D < 0.50 A, and between 0.50 
and 1.44 A; and again for 17 O-type stars. These computations 
with a few others for comparison are summarized in Table 5. 

What is the proper interpretation of these results, which, taken at 
face value, indicate that \ 6284 has a closer connection with space 
reddening than the D lines? Is a physical relationship between space 
reddening and intensity of \ 6284 to be inferred? Since the D lines 
arise from absorption by gaseous sodium atoms, while space redden- 
ing is probably caused by finely divided dust, a close physical de- 
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pendence in this instance is unlikely. Were both sodium gas and the 
dust uniformly distributed through space, a very high correlation of 
I'D and color excess would nevertheless be expected. The moderate 
values of the coefficients recorded in Table 5 indicate considerable 
irregularities in distribution. May we conclude from the stronger 
correlation between color and J’ \ 6284 that the source of this line is 
in some way connected with the dust particles? Errors of measure- 
ment are larger for \ 6284 than for the D lines and would therefore 
tend to weaken the correlation relative to that of the D lines rather 
than to strengthen it. On the other hand, the approach to satura- 
tion may give more scatter to the values of I'D, the effect being 


TABLE 5* 
CORRELATION COEFFICIENTS 

Objects No. E,I’D E, I'd 6284 I'D, I'd 6284 

Stars having /All...... 85 | +o.61+0.05 | +0.81+0.03 | +0.68+0.04 
complete data Weak D.. 45 .40+ .08 {S47 507 .46+ .08 
Strong D. 40 81+ .04 .46+ .08 

17 | +0.44+0.13 | +0.61+0.10 | +0.39+0.14 


* FE = Color excess (Stebbins and Huffer, and Becker); J’D = intensity 0.5(D2+ D1); I’ 6284 = 
“smoothed” intensity of \ 6284. 


negligible for \ 6284. Even for the D lines, however, the effect can- 
not be very great, unless possibly for extremely intense lines, be- 
cause intensity is, on the average, nearly proportional to distance."® 


ORIGIN OF THE LINES 


Because of the very long intervals between collisions resulting 
from the low density, particles in interstellar space capable of emit- 
ting light must spend practically all their time in the quantum level 
of least energy. Even levels ordinarily called ‘‘metastable’’ prob- 
ably do not persist long enough to cause any important exceptions 
to this rule. Interstellar absorption lines should therefore be reso- 
nance lines. Atomic lines of this kind, with the possible exception 
of those belonging to very rare atoms, should be well known in the 
laboratory. The difficulty of identification of the new interstellar 


19 The mean distance-intensity relationship is X See Mt. W. 
Contr., No. 569; Ap. J., 86, 28, 1937. 
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lines in the visual region may thus be taken as an argument against 
an atomic origin. 

The width of the lines, moreover, seems incompatible with an 
atomic origin. Since the lower level of an interstellar absorption line 
must have an extremely long life, the line would be very sharp. The 
upper level might conceivably be one of very short life, perhaps with 
energy in excess of ionization, but this situation seems improbable 
as it would call for a lower level of high energy with incredibly great 
metastability. 

An approximation to a (hypothetical) damping constant for the 
line \ 6284 may be obtained as follows. Since the line is relatively 
shallow even for a considerable total intensity, the equivalent width 
in angstroms may be written 

mc? 


W = 10° Nf. 


For an equivalent width of 1 A, this equation gives Vf = 2.9 X 10”. 
The optical thickness is 
=Nfo, 


where a, the atomic absorption coefficient, may be written 


6 


6 being the quantum-mechanical analogue of the classical damping 
constant. At the center of the line 

ro me 

Table 3 shows that for an equivalent width of 1 A, \ 6284 has a 
central absorption of 20 per cent. On account of the width of the 
line, this value is probably not seriously affected by photographic 
effects or by the finite purity of the spectrum. Hence we may write 


= 0.206, 


from which 


2.423 % 10” 
T = 0.22 = 
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Hence 6 = 1.1 X 10", which is about two thousand times the value 
of the classical damping constant 

= — = 7 
ae 5.7 X 107 (for \ 6284) . 
If \ 6284 is produced by an atom whose thermal motions are 

responsible for the width, the atomic mass must be very small. If 
the optical thickness at the center of the line is 0.22, that at distance 
Ad from the center may be taken to be 


Tar = 0. 


Our measurements show that, for an absorption of 5 per cent, 
Ad = 3.4A. Hence = 2.8 A. Since 


it follows that 
m=1.5X 10% XT. 


For the lightest atom known, that of hydrogen, T = 1,100,000°, an 
inadmissibly high value, while for heavier atoms T is still greater. 
For T = 40,000°, which is probably an upper limit,?° m = 6 X 10°”, 
which is only one-thirtieth the mass of the hydrogen atom, while 
for T = 10,000° it is 1.5 X 10°, or only about twenty times the 
mass of the electron. These considerations make it very improbable 
that \ 6284 is an atomic line broadened by Doppler effect. 

A second hypothesis is that the lines are actually bands or por- 
tions of bands of molecular origin. Under the prevailing excitation, 
molecules would be concentrated in a few of the lower quantum 
levels, and the band structure might thus differ from that observed 
in the laboratory. The low density of interstellar space favors both 
ionization and dissociation, but since neutral atoms (of sodium) are 
known to be fairly abundant, one can scarcely decide a priori against 
the existence of molecules. Among the writers who favor the mo- 
lecular hypothesis, two have made specific suggestions. P. Swings” 
has called attention to the fact that \ 6613.9 and \ 6283.91 fall 
approximately in the calculated positions of bands of carbon di- 
oxide. The second coincidence is not very close (AA = 3.0 A), but 


20 Mt. W. Contr., No. 570; Ap. J., 86, 44, 1937. 7! M.N.*, 97, 212, 1937. 
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Swings thinks it may be within the uncertainties of measurement 
and calculation. A spectrogram of x? Orionis taken by Dr. Theodore 
Dunham, Jr.,”” which includes the region of the ‘‘Venus” bands at 
d 7820 and A 7883 and fails to show definite absorption at or near 
these wave lengths, gives a negative answer to the test suggested 
by Swings. 

M. N. Saha? has recently suggested that \ 6284 ‘‘appears de- 
cidedly to be due to molecular sodium, Na,, due to the transition 
(Na,: — B’S, n” = 0, n’ = 8)” and that 25780 
“may be provisionally identified with a line of NaK (n”’ =o, 
n’ = 5).”’ Detailed calculations of the structure of these bands and 
of the relative intensities of neighboring related bands will probably 
be required before the identifications can be conclusively established. 

A third possibility is that the unidentified lines are produced by 
small solid particles. This hypothesis is suggested by the fact that 
the correlation between color excess and the intensity of \ 6284 
seems to point toward a physical relationship between space red- 
dening and the formation of the spectral line. At room temperature 
only a few solids are known to produce narrow absorption lines, 
but at very low temperatures many, although not all, substances 
probably have narrow lines which may be thought of as displaced 
atomic lines.*4 It is therefore conceivable that at temperatures a 
few degrees above absolute zero, the interstellar dust which causes 
the reddening of star light might also produce certain detached 
lines. This surmise, if correct, offers the future possibility of a 
partial chemical analysis of dark nebulae. A theoretical investiga- 
tion of the formation of absorption lines by solids in the form of 
finely divided dust would be valuable, while data on lines produced 
in the laboratory by solids at low temperatures should be examined 
for possible coincidences in wave length. 

An alternative interpretation of the data might be that the lines 
have their origin in a gas which is in some way dependent on the 
dust particles and thus tends to occupy the same regions of space. 

CARNEGIE INSTITUTION OF WASHINGTON 


Mount WILSON OBSERVATORY 
August 1937 


22 Who kindly permits reference to it here. 3 Nature, 139, 840, 1937. 
24 See, e.g., Spedding and Nutting, J. Chem. Physics, 2, 421, 1934. 
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ROBERT B. KING AND ARTHUR S. KING 


ABSTRACT 


The relative f-values for 115 lines in 19 multiplets of Fe 1 between \ 2912 and d 5455 
and for 227 lines in 66 multiplets of Ti 1 between \ 3123 and A 6743 have been derived 
from measurements on electric-furnace absorption spectrograms. The method is es- 
sentially the same as that described in Mt. Wilson Contr., No. 528, which gave results for 
75 lines in 9 multiplets of Fe 1. 


All the lines arise from the lowest terms in the atoms, ranging in excitation potential 
from o.0 to 1.5 volts for Fe 1 and from 0.0 to 2.5 volts for 771. 

The relative f-values of lines in multiplets of Fe 1 and 771 have 
been measured with the electric furnace as the source of the ab- 
sorption-spectrum excitation. The method used has been described 
in detail in an earlier paper on relative f-values for lines of Fe 1.? 
This paper will be referred to in the following as “‘(1).”” The method 
will be summarized briefly here, and modifications made in the 
present work will be pointed out. 

The elements have been vaporized in the furnace at known tem- 
peratures, and their spectra observed in absorption. The relative 
total absorptions of the lines appearing on the linear (faint-line) por- 
tion of the curve of growth were measured and, with the aid of this 
curve, converted into relative intensities for the given temperature. 
The relative f-values are then obtained by the application of Boltz- 
mann factors to the observed relative intensities of lines arising 
from different energy-levels. Boltzmann factors may be applied be- 
cause of the existence of a close approximation to thermal equi- 
librium at a known temperature in the furnace. This condition was 
tested in (1) and found to hold within the limits of experimental 


error. 
OBSERVATIONS AND PHOTOMETRY 


The absorption spectrum to be studied was obtained by passing 
the light from a 50-amp. Philips cinema lamp with quartz bulb 
* Contributions from the Mount Wilson Observatory, Carnegie Institution of Washing- 


ton, No. 581. 
"R. B. King and A. S. King, Mt. W. Contr., No. 528; Ap. J., 82, 377, 1935. 
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through the furnace tube and to the slit of the 15-foot concave- 
grating spectrograph. The newer hood furnace was used in place 
of the cylindrical furnace employed in the earlier work because a con- 
stant temperature was more readily maintained. 

The spectrum to the violet of about \ 5400 was, as a rule, photo- 
graphed in the second order, while the first order was used to the 
red of this limit. Eastman IV-O and IV-F plates were used, the 
fine-grained, high-contrast emulsion being most satisfactory for the 
photometry of faint absorption lines. 

In the previous work all the exposures were made at a furnace 
temperature of 2100° C, the intensities of the lines being altered for 
different exposures by varying the charge of material in the furnace. 
This alteration of intensity is necessary in order that a large in- 
tensity range may be treated, the linear portion of the curve of 
growth covering only a short range. The change in the strength of 
the lines can, however, be more readily controlled by changing the 
temperature, as well as the charge of material, in the furnace for 
different exposures; in the present work, temperatures ranging from 
1800° C to 2600°C were used. The low temperatures bring the 
strongest lines to the linear portion of the curve of growth, and the 
high temperatures the fainter lines. Practically all the Fe and Ti 
lines of temperature classes 1 and 11 and many of class 11 were thus 
brought onto the straight portion of the curve. 

Temperature measurements were made on the inner walls of the 
furnace tube by means of an optical pyrometer, checked for accuracy 
by means of a tungsten ribbon lamp calibrated by the Bureau of 
Standards. Correction was made for absorption by the furnace 
window. The uncertainty in the temperature readings was usually 
less than 20°C. The actual uncertainty, of course, also depends 
upon the uncertainty of the temperature scale at these temperatures. 

Calibration marks were impressed on the plates by means of a 
ribbon-filament lamp and a step-slit, as described in (1). Plates were 
developed in X-ray developer at 18° C, as before. 

Tracings of both furnace and calibration spectra were made with 
the photoelectric microphotometer designed by Dunham. The high 
resolving-power and sensitivity of this instrument are well suited to 
photometry of faint absorption lines. 
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REDUCTION OF OBSERVATIONS 

Characteristic curves constructed for each plate at wave-length 
intervals of about 100 A were used to derive the equivalent widths 
of the absorption lines in angstrom units. In the work reported in 
(1), the areas under the density contours of the lines were converted 
into intensity areas, which were then measured. For the present 
work a shorter method has been adopted. It happens that the 
Doppler width of lines of Fe and Ti at the upper limit of the linear 
portion of the curve of growth for the temperatures used is less than 
the width of the slit image formed on the plate even in the second 
order. That is, with a slit-width of 0.015 mm the faint lines falling 
on the linear portion of the curve of growth appear on the micro- 
photometer tracings as small triangles with the same base-width. 
The altitudes of the triangles are therefore proportional to the 
equivalent widths of the lines. The base-width varies directly with 
the wave length and, in angstrom units, is very nearly twice as large 
in the first order as in the second, owing to the increased resolving- 
power of the second order. Mean values of the base-width as a func- 
tion of the wave length were determined for the first and second 
orders and used together with the measured altitudes of the triangu- 
lar contours to obtain equivalent widths. 

The approximate linearity under our experimental conditions of 
the ‘“‘straight’”’ portion of the curve of growth for faint lines was 
verified in (1). In the present work an empirical curve of growth 
was used to reduce the equivalent widths to relative intensities. 
This curve was constructed by plotting the equivalent widths, ob- 
tained by treating the line contours as triangles, against Harrison’s? 
arc-intensity measurements for several multiplets of 771 as described 
in (1). The resulting curve differed only very slightly from the 
theoretical linearity. 

The total absorption A; of a line formed by passing a continuous 
spectrum through an absorbing column of gas has been defined? as 
the ratio of the energy absorbed by the line to the incident energy: 


2G. R. Harrison, J. Opt. Soc. Amer., 17, 389, 1928. 
3 R. Ladenburg and F. Reiche, Amn. d. Phys., 42, 181, 1913. 
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On the linear portion of the curve of growth the total absorption A; 
in a frequency interval dv is directly proportional to the number of 
atoms absorbing the line: 


« Nf, 


where N is the number of atoms in the lower energy-level involved 
in the transition and f is the f-value, or “oscillator strength” of the 
line. Then 


Hence, if d\ is measured in angstrom units, the quantity 47 dd is 
the equivalent width EW, the fraction of an angstrom unit of the 
neighboring continuous spectrum totally absorbed by the line. So 
we can write 


« Nf. 

We are concerned here only with the relative numbers of atoms in 
the various low levels. These depend upon the temperature of the 
gas, the distance of the levels above the normal state, and the sta- 
tistical weights, g, of the levels. If the gas is in thermal equilibrium, 
a Boltzmann factor may be used to allow for the temperature dis- 
tribution of the atoms among the low energy-levels. Hence, 

he 
Nae 


where Av is the wave-number difference between the low level from 
which the absorption line arises and the lowest energy-level in the 
atom, and T is the absolute temperature. We have then 


and 


| 
Are Nf. 
\ 
he 
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EW 
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The procedure followed was to derive from the linear portion of 
the empirical curve of growth numbers corresponding to the relative 
intensities of the lines on each spectrogram. Then, since a variety of 
temperatures was used, Boltzmann factors were applied. Intensities 
from different spectrograms were reduced to the same scale by means 
of factors obtained from lines appearing on two or more spectro- 
grams, and all the data for each line averaged to obtain the adopted 
mean values.‘ These intensity numbers, when divided by }’, are rela- 
tive gf-values. It has been pointed out that to correspond to true 
relative f-values these must be divided by the statistical weight, 
g = (27 +1), of the lower level involved in the transition.’ Since 
this was not done in (1), the values there published are really gf- 
values. Since relative gf-values are also convenient for astrophysical 
applications, they, as well as the relative f-values, are listed in the 
tables of results. 

RESULTS 

Fe 1.—Relative gf- and f-values, the former on an arbitrary scale, 
for 115 lines in 19 multiplets of Fe 1 between \ 2912 and A 5455 are 
given in Table 1. For completeness the results for the 75 lines in 9 
multiplets between A 3649 and A 4325 published in (1) are also in- 
cluded in the table. The number of spectrograms used in the means 
(last column of table) indicates roughly the weight of the results. 
Values which are particularly uncertain because the lines were very 
faint on the spectrograms are marked with an asterisk. A few of the 
strongest lines whose measures were regarded as uncertain are also 
thus indicated. 

All lines measured arise from the a‘D, a‘F, and aF levels, the 
three lowest even levels in Fe 1, whose excitation potentials are 0.0, 
1.0, and 1.5 volts, respectively. Only lines arising from these three 
low terms are observed in absorption at the temperatures used. 

Of the 29 multiplets listed, only 7 are completed. In most cases 
the fainter members were not observed. For multiplets Nos. 172 and 
203, however, the strongest lines were too strong to appear on the 

4 The method follpwed in (1) differed from this procedure in that, since all measures 


were made at a single temperature (2100° C), a single set of Boltzmann factors for this 
temperature was applied after the measures were reduced to the same scale. 


5D. H. Menzel and L. Goldberg, Ap. J., 85, 40, 1937. 
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TABLE 1 
RELATIVE gf- AND f-VALUES FOR LINES OF Fe I 
Mult. No. r Multiplet af f No. of Plates 
2912.16 10 2 
2929 .O1 a5D;—y5F3 13.5 1.9 2 
2930.90 a3D,—y5F4 25 2.8 6 
2941 .35 aSD.— y5F? 12 2.3 I 
2947 .88 as3D,—y5F3 28 4.0 6 
2053-94 a5D,— y5F3 28 5.6 6 
2957.37 a5D,—y5F? 22 7.4 6 
2965.26 aS5Do— y5F? 17 17 5 
2966 .go asD,—y5Fs 43 4.7 7 
2970.10 aSD,—y5F3 27 9.0 6 
2981.45 a5D;,—23P3 14 2.0 8 
2986.46 aS3D,—23P? 1.2 5 
3007 . 29 II 6 
3024.04 a3D,—23P2 16 5.2 8 
2983.57 aSD,—y5D; 36 4.0 6 
2904.44 asD;—ysD2 43 6.1 7 
3001 .00 46 Q.2 6 
3008.14 35 12 6 
3017 .63 16 5.2 6 
3020.49 asSD,—ysD? 29 5.8 6 
3020.64 51 5 
3021.07 aSD,—ysD3 63 9.0 6 
3025.85 asDo—ysDi 34 34 6 
3037.39 asD,—ysD3 45 15 6 
3047.61 a5D,— 45 9.0 
3059 .09 38 5-4 7 
2987 . 30 asF,—x5F3 180 20 5 
2999 52 340 31 5 
3003 .04 a5F,—x5F3 190 27 5 
3009.58 320 36 5 
3016.20 asF,—x5FY 160 32 5 
3018.99 a5F,—x5F$ 310 44 4 
3026.47 asF,—x5F2 270 54 5 
3031.64 aSF,—x5FY 220 73 5 
3040. 43 asF,—x5F3 160 18 5 
3041.75 ask, 230 33 5 
3042.03 aSF,;—x5F3 130 44 4 
3042.67 aSF,—x5F3 220 44 5 
3057.45 500 46 5 
3067 . 25 aSF,—x5D§ 520 58 5 
3075.73 asF;—xsD2 490 7° 5 
3083.75 aSF,—x5Dz 380 76 5 
3091.58 asF,—x5Do 270 go 5 
3099.90 asF,;—xSDi 370 123 5 
3099.97 asF,—xsDj 310 34 5 
3100. 31 aSF,—x5D3 390 78° 4 
3100.67 aSF;—x5D§ 350 50 4 
3116.64 asF,—x5D3 95 32 3 
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TABLE 1—Continued 
Mult. No. r Multiplet af - No. of Plates 
g1—Cont.| 3125.66 aSF,—x5D¢ 120 24 3 
2534.11 asF;—xsDj II 3 
3143.24 0.69 0.077 3 
3180.76 | 0.68 4 
3184.90 asD;—2F3 6.1 0.87 5 
3193.22 asD,—2F§ 6.9 0:77 5 
3199.53 aD, —2F3 2.9 0.97 5 
3214.40 1.10 5 
3230. 23 asD;—23F 6.1 0.87 5 
3191 .66 asD,—2D3§ 0.57 5 
3196.94 as5D,—2D? 4.1 0.59 4 
3200.79 asD,—2D? 0.86 0.17 4 
3229.13 aD,—2D? 2.3 5 
3234.62 32.8 0.54 4 
3245.98 a5D,— 23D? 1.4 5 
3205.06 aSD,—23D$ 0.66 5 
3452.28 aSF;—y3F5 66 9.4 2 
3476.71 a5Do—z5P? 65 65 I 
3407 .84 23 I 
3520.07 a’D.—z5P$§ 12 2.4 I 
3513.82 asF;—z3G§ 260 24 6 
3521.26 asF,—2G3 310 35 6 
3520.17 aSF;—23G3 180 24 5 
3558.52 aSF,—23G3 660 132 6 
3565.38 1800 250 <2 
3570.10 aSF,—z3G¢ 2600 280 I 
3585.32 480 69 3 
3585.71 260 29 4 
3586.99 aSF,—z5G¢ 550 110 6 
3589.11 5° 4.5 3 
3608 . 86 aSF,—25G$ 2000 660 2 
3618 ..39 aSF,—z5G3 2900 590 2 
ale 3644-46 asF;—25G§ 2000 270 I 
3647.85 asF,—z25G5 1600 170 I 
3679.92 aD,—2aF3 61 6.8 16 
3683.05 aSD,—25F2 9.0 1.3 9 
3705.57 a’D,—25F3 IIo 16 13 
3707 .83 as8D,—25F 13 2.6 9 
3719.94 asD,—25FS 530* 59* 7 
3722.56 a’D,—2ZF? 120 24 13 
3733-32 aSD,—25F? g2 31 14 
3737-13 asD;—25F 450* 64 8 
3745.50 asD,—25F3 340 68 8 
3745.90 110 110 16 
3748.26 asD, 190 63 II 
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Mult. No ny Multiplet gf f No. of Plates 
3687.46 asF,—y5F4 430 390 7 
3709. 25 asF 710 79 10 
3727.62 aSF;—y5F3 760 110 9 
3734.87 asF;—y5FS 4700 43 12 
3743 -34 a5F,—y5F? 560 110 8 
3749-49 ask, —ysFy 3700 410 13 
3758.24 a5F—y5FS§ 2800 400 13 
3763.79 a5F,— y5F3 1800 360 II 
3767.19 asF,— 1300 430 10 
3787.88 — y5F3 500 170 4 
3795.00 590 120 7 
3798.51 a5F,— y5F 260 29 5 
3799-55 asF;—ysF 500 72 6 
BOBS 3786.68 12" 7 
3790.10 asF,—23P¢ 04 19 7 
3812.97 asF,—23P3 330 47 8 
3814.53 24* 8.0* 8 
3850.82 aSF,—z23P3 100 20 8 
308.......} a3F,—y3D§ 5200 580 8 
3827.83 aF;—y3 D2 4700 670 8 
3841.05 aF,—y3D? 3600 720 8 
3888 52 a3F,—y3D? 1200 240 
3902.95 aF,—y3D§3 1800 260 8 
3820.43 ySD§ 3100 280 16 
3825.88 aSF,—y5D§ 2200 250 16 
3834.23 a5F;—ysD$ 1500 210 II 
3840.44 aSF,—ysD? 870 170 II 
3849.97 asF,—ysDo 510 170 9 
3865 .53 aSF,—ysD? 420 140 9 
3872.50 a5F,—ysD3 410 82 II 
3878.02 420 60 II 
3887.05 asF,—ysD3 270 30 10 
3898.01 aSF,—ysD3 53 19 5. 
3917.19 asF,—ysD§ 42 8.4 4 
3940.88 16* 2:3" 5 
3824.44 asD,—25D$ 87 9.7 12 
3850.37 a5D;—z5D? 16 12 
3859.91 asD,—25D§ 350 39 6 
3878.56 as’D,—z5D? 88 17.6 12 
3886. 28 160 23 10 
3895 .66 aS5D,—z5D¢ 52 17 13 
3899.71 asD,—z5D$ 70 14 13 
3906. 48 as’D,—z5D? 19 6.3 10 
3920.26 asD.—z5D? 48 48 13 
3922.91 58 $33 15 
3927.92 a5D,—25D3 70 23 16 
3930.30 a5D,—25D$ 78 16 16 
3969. 26 aF,—y3F$3 1500 170 12 
4005.25 aF;—y3F2 1200 170 II 
4045.82 6000 670 8 
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TABLE 1—Continued 
Mult. No. as Multiplet af f No. of Plates 
353—Cont.| 4063.60 aF;—y3F3 4300 610 13 
4071.75 aF,—y3F? 3800 760 12 
4132.06 aF,—y3F3 1100 220 9 
4143,87 @F;—y3F3 1400 200 II 
WES 4202.03 890 99 18 
4250.79 ak 880 130 16 
4271.76 a@F,—2G5 2450 270 19 
4307.91 3300 470 16 
4325.77 a3F,—23G$ 3800 760 18 
4206.70 0.60 0.086 ir 
4210.19 aSD,—z'P§ 2.14 0.24 16 
4258.32 aSD,—z'P$ 0. 38* 0.076* 4 
4291.47 0.45 0.064 12 
4204.13 ak 380 42 6 
4337-05 a3F,—z5G¥ 130 19 2 
4383 .55 —25G5 5800 650 5 
4404.75 asF,—25G5 3000 430 10 
4415.13 a3F,—25G3 1000 200 7 
4375-93 4.4 °.49 4 
4389. 25 0.24* 0.034 3 
4427.31 4.3 0.61 10 
4435.15 aSD,—z'F? 0. 24* 0.048* 2 
4461 .66 2.9 0.58 12 
4466.56 0.31 0.10 3 
4482.18 1.6 0.54 12 
4489.74 0.63 0.63 7 
BOS. 4602.95 a3F,—ysFS 120 13 I 
5012.07 asF;—z5F5 II 1.0 3 
110.41 asD,—z’D§ 0.76 0.084 7 
5166.29 a5D,—27D§ 0.35 0.039 4 
5168.90 0.53 0.076 5 
5204.59 0. 26* 0.052* 3 
5167.49 aF,—z3D§ 260* 29* 5 
5270.36 180* 30* 
5269.54 aSF;—zsD} 170 15 9 
5328.04 asF,—25D§ 120 13 9 
5371.50 aSF;—z5D$ 80 II 9 
5397-13 asF,—z5D3 41 4.6 6 
5405.78 aSF,—z25D¢ 58 12 7 
5429.70 asF;—z5D$ 51 7-3 7 
5434.53 34 II 6 
5440.92 aSF,—z5D2 43 8.6 8 
5455.62 aSF, —25D¢ 35 12 
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linear portion of the curve of growth on any of the spectrograms. 
Both these multiplets arise from the a‘D term, and their gf-values 
are probably somewhat larger than those for the strong lines of 
multiplet No. 257. New measures have been included in the results 
for Nos. 257, 353, 418, and 412, which were published in (1). The 
relative gf-values for a few of the lines of these multiplets have been 
_ slightly revised by means of the additional data. These lines are 
AA 4045, 4063, and 4071 in No. 353; AA 4202, 4250, 4271, and 4307 
in No. 418; and AX 4216, 4258, and 4291 in No. 412. 

All the Fe lines listed were observed in the second order, the first 
order being used only to connect the lines of the \ 4400 and X 5200 
regions in reducing them to the same scale. 

Tit.—Results for 227 lines in 66 multiplets of 7z1 between 
d 3123 and A 6743 are listed in Table 2. The relative gf-values are 
on an arbitrary scale, which has, of course, no relation to that for 
lines of Fe 1. 

As in the case of Fe 1, all the observed lines arise from the lowest 
terms, namely, a3P, b3F, a5P, a3G, z5G°, b'G; 
and a'P, ranging in excitation potential from 0.0 to 2.5 volts. 

In only 33 of the 66 multiplets are the observed intensities com- 
plete. The missing lines are all weaker members which did not ap- 
pear in absorption at the temperatures used. 

With the exception of No. 112, which includes both first- and 
second-order measures, all the multiplets to the red of \ 5100 were 
measured in the first order. Because of the greater dispersion and 
resolving power in the second order, measures made there have less 
uncertainty than those made in the first order. 


SOURCES OF ERROR 


Sources of error were discussed in (1). It may be emphasized 
here, however, that in addition to accidental errors inherent in the 
photographic photometry of faint absorption lines and possible 
systematic errors due to temperature measurements (and uncer- 
tainty in the temperature scale at high temperatures), there are 
opportunities for the introduction of systematic errors in reducing 
the measurements from different spectrograms to the same scale. 
Wide ranges in intensity and in wave length were covered, and in 
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TABLE 2 
RELATIVE gf- AND f-VALUES FOR LINES OF 77 I 
Mult. No. r Multiplet gf f No. of Plates 

BO Ave 3123.07 a'D,—w'P? 560 112 5 
eee 3141.51 a'D.—x'P? 440 88 5 
3186.45 1200 240 3 
3192.00 aF;—w3G4 1750 250 3 
3199.91 a3F,—wiG? 2300 256 2 
3203.82 aF;—wsG3 123 17.6 4 
3214.23 a@F,—w3G3 134 15 4 
3226.22 ‘3° 0. 48* 5 
3204.87 aP,—s3D$ 280* 93* 5 
3205.85 a@F,—v3D? II 2.2 5 
3221.14 a@F;—v3D$ 4.0* 6.57" 
3222.74 aF;—v3D$ 12.2 5 
3243.80 a@F,—v3D$3 10.5 
3226.11 2>G8—g5Ge 16000* 1200* I 
3292.00 a'D.—x'F$3 1060 210 5 
3209.41 a'D,—v3G$ 240* 48* 5 
3308. 38 aP,—t3D? 560 560 5 
3309.49 asP;—t8D3 1300 440 5 
3314.42 a3P,—t3D§} 2200 440 4 
3314.50 a3P,;—t3D? 520 170 3 
3321.58 a3P,—t3D$ 440* 88* 5 
3341 .87 aF,*x3G3 1080 220 5 
3354.04 1570 220 5 
3301.00 85 120 9 
3371.44 aF,—x3Ge 2050 230 5 
3379.20 142 16 5 
3385 .66 aF,—x3G$ 84 9.3 9 
3358.26 86 17.2 8 
3361.27 @F,—w3D§3 120 172 5 
3370.42 a@F,—wsD? 490 98 q 
3377-57 aF,;—wD3 690 99 7 
3385.93 775 86 7 
3342.70 a3F,—x5D$ 6.4 5 
3348.52 a3F,—x5D¢ 7.4 Es 5 
3352.92 aF;—x5D§ 29 4.1 
3301.82 a3F;—x5D3 48 6.9 5 
3377.48 @F,—xsD§ 570 63 4 
3382.30 asP,—wD§ 660 130 5 
3390.67 aP,;—wD$ 460* 150* 5 
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Mult. No r Multiplet gf f No. of Plates 
BOM 3480.53 550 110 5 
3493.27 6:9 0.96 5 

3506.64 @Fy—ysF§ 22 2.45 5 

3603 . 86 a3F,—25P3 0. 48* 5 
3006. 81 6.2 0.89 II 

3626.09 6.x 0.87 II 

3635.20 @F,—Z2P3 26 2.9 II 

3604. 30 9.1 II 
3010.16 a'D,—y'P? 850 170 II 
3635.47 aF,—y3G$ 1600 320 7 
3042.68 1900 270 7 
3653.49 aF,—y3G5 2500 280 6 

3658.14 117 17 10 

3671.66 a@F,— 121 13.5 10 

3687.38 @F,—y3G§ 9.6 II 

3637.97 aF,—xD§ 19 3.8 II 
3046.19 aF,—x3D$3 40 8.0 14 

3654.58 62 12.4 17 

3660.62 aF;—x3D¢ 58 8.3 16 

3668.95 aF;—x3D$ 71 10.1 17 

3689 . 89 aF,—x3D3 78 8.7 17 

3702.29 a3Po— 250* 250* 3 
3709.95 610 200 10 

3725.12 a3P,—y3S? 970 194 II 

3717.39 81 16 9 
3722.58 @F,;—xF4 95 14 5 

7 3729.77 690 140 8 
3741. 1000 150 4 

3752.87 BF 1700 190 4 

3753-63 122 17:5 10 

3771.64 @F,—xXF§ 129 14.3 10 

Se 3724.59 a'G,—x'G§ 3400 380 8 
3786.03 a'D,—z'P? 1450 290 6 
3900.96 40 6.6 7 

3914.33 vF,—ysD; 75 8.3 9 

3934.23 aF,—ysD3 9.6 10.7 6 

3898. 49 a@F,—y3D3 11.6 6 
3924.51 aF,—y3D§3 180 26 6 

3929.87 aF,—y3D$ 132 26 6 

3948.66 @F,—y3D? 590 120 5 

3956.28 aF;—y3D3 610 87 5 

3958.21 1120 125 5 
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TABLE 2—Continued 
Mult. No. r Multiplet gf - No. of Plates 
3904.77 a'D,—y'F$ 2300 400 6 
3914.72 4.6* 0.92* 6 
3921.42 aF,—23P?2 44 8.8 6 
3947-75 167 24 6 
3902.85 105 21 6 
3064.27 105 15 6 
12 3981.77 770 154 7 
3989.77 1090 156 7 
39908 .65 1700 190 7 
4008 .94 151 22 8 
4024.57 @F,—y3sF3 164 18 6 
3982.54 78 15.6 II 
4009 .68 ak; —z5S$ 32 4.6 8 
4055.02 290* 290* 6 
4060.27 370" 120* 6 
4004.22 a3P,;—x3P¢ 250* 83* 4 
4065.11 aP,—x3P> 290* 97* 6 
4078.47 1070 210 4 
4082.46 a3P,—x3P? 320 64 ig 
4076.38 a3F;—z'G§{ 2.0* 0. 29* I 
4112.73 40 45 6 
209...% 4272.44 asF;—x5D§ 70* 10* 2 
4274.60 asF,—x5D§ 106* 21* 5 
4281.40 aSF,;—x5D? 110 37 
4286.01 aSF;—x5D3 780 110 10 
4287.41 asF,—x5D{ 780 86 10 
4289.08 aSF,—x5D3 950 190 10 
4290.93 aSFi—xSD? 680 230 10 
4205.75 asF,—x5Do 740 250 10 
4298.67 1310 260 10 
4300.55 adF;—x5D3 2100 300 4 
4301.08 asF,—x5D§ 2600 290 4 
4305.91 46c0 420 3 
4299.64 asF;— 240 34 8 
asF,— wD? 
4326.35 asF;—w3D3 140* 20* 3 
4465.80 aSP,—y5P$ 1400* 280* 2 
4471.24 asP, — y5P$ 1500* 500* 2 
4481.27 a5P;—ysP$ 2800* 400* 3 
4489.10 1800* 360* 2 
172..44..| 4512.74 asF,—ysF2 550 61 13 
4518.03 750 107 12 
4522.80 asF,—ysF$ 740 148 12 
4527.32 asF,—y5F? 520 173 13 
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Mult. No. r Multiplet gf f No. of Plates 
172—Cont.| 4533.25 a5F;—ysF§ 5100 460 6 
4534.78 asF,—ysF 3500 390 5 
4535.58 1900 270 4 
4535.92 aSF,— y5F3 1180 240 10 
4530.00 asF,—y5F¢ gIo 300 12 
4544.70 a5F,—y5FY 500 100 13 
4548.72 asF,—ysF3 735 105 14 
4552.45 755 84 12 
4555-50 550 50 13 
4562.63 aF;—z'D? 3.6 0.51 4 
4617.27 a5P;—wsD§ 3400 490 5 
4623.10 asP,—w5D§ 1700* 350* 5 
4639.36 a5P,—wsD? 1500* 300* 2 
4639.66 a5P;—wsD§ 1300* 180* I 
4650.47 @F,—23G3 70 14 8 
4667.50 gI 13 8 
4681.91 aF,—23G§ 132 14.7 8 
4093.68 @F;—2G3 2.6 0.37 5 
4715.31 aF,—23G4 2.4 0.27 5 
4675.13 a3P,—x5D§ 14* 2 
res. J. 4091.34 asP,—wsD§ 210 42 5 
4698.79 aP,—wsD? 160 54 4 
4710.19 a3P,— wD? 80* 80* 4 
4722.63 a3P,—wsD? 66* 22* 2 
4840.88 a'D,—y'D? 460 92 6 
338..1$ <2). 4885.09 aG;— 3400* 310* 4 
4899 .93 asG,—yH§ 2600* 290* 3 
4913.63 2500* 360* 2 
4981.75 y5G8 5200 470 6 
4991 .08 aSF,— y5G§ 4000 440 6 
4999.51 asF;—ysG4 3150 450 5 
5007.22 a5F,—ysG$ 2100 420 4 
5014.28 aSF,—y5G$ 1600 540 4 
5016.17 aSF;— 460 42 6 
5020.04 asF,—y5G; 640 71 6 
5022.87 aSF;— y5G$ 590 84 6 
5024.85 asF,—y5G$ 410 82 6 
4997.10 aF,—z23D? 12 2.4 6 
5009.65 aF;—2D§ 8.9 6 
5014.19 94 19 6 
5039.96 @F,—2D3 135 19 6 
5064.66 180 20 5 
2500 280 5 
5030.47 1950 280 5 
5038.41 bF,—wiG$ 1550 310 5 
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TABLE 2—Continued 
Mult. No ny Multiplet gf if No. of Plates 
5113.45 b3F;—v3D$ 420 60 4 
5145.47 bsF,—v3D; 670 74 10 
5147.48 BF, 2.9 12 
5152.18 14.4 12 
5173-74 aF,—2F> 135 27 5 
5192.98 @F;—2F3 180 20 5 
5210.39 BF, 240 25 4 
5219.72 | II 
6262. @F,—2F3 4.8 0.53 5 
5206.09 a'P,—w'P? 16500 5500 I 
5238.58 asF;—ysD§ ax 8 
5240.56 aF,—ysD3 1.95 7 
5250.94 a5F;—ysD3 5 
26). 5205.07 1100 100 3 
5283.45 a3G,—v3F$ 930 100 2 
5207.29 aG;—v3F3 860 120 3 
5282.38 a3P,—xD? 42 14 
5284.39 30 30 2 
5295.78 asP,—x3D3 68 13.6 10 
98...25.. 5306.63 asF,—y3F$ 11* I 
5389.15 asF,—y3F2 9. 2 
5408 .93 0. 38* 0.076* 2 
5420.27 2255 0.36 
5440.67 aF;—z2D? 0.144 
04. 5474.28 b3F,—x3G¢ 140* 16* 4 
88. ./0G.. 5481.90 b3F,—w3D$ 140* 28* 4 
5512.54 b3F,—w3D5 840 93 10 
5514.36 bF,—w3D¢ 640 128 
5514.54 bF,—w3D$ 860 123 
92.../9.7. 5490.16 bF,—x5D§ 300 33 7 
5644.14 b'G,—y'G§ 3300 370 I 
5866.45 265 53 9 
5899. 29 aP,—y3D? 157 52 IO 
5922.11 aPo—ysDi 75 75 9 
5041.73 asP;—ysDi 57 19 8 
5880.30 30 10 I 
5918.56 50 7 
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TABLE 2—Continued 


Mult. No. ny Multiplet gf f No. of Plates 
5953.16 aG;—2H¢ 1100 100 3 
5965.83 a3G,—2H8 1000 110 4 
5978.54 aG,—23H4 850 120 2 
6085 . 22 a3P,— 238% 64 21 10 
6126.21 a3P,—28¢ 82 16 10 
6258.10 bF;—y3G} 670 96 10 
6258.70 b3F,—y3G¢ 875 07 10 
6261.10 bD3F,—y3G3 560 EE2 10 
6546.26 120* 24* 2 
6554.23 b3F;—x3F3 120* 17” 4 
6556.08 155 17 4 
6599.12 a'D,—2'F$ 15 4 
ee 6743-14 a'D,—z'D? 48 9.6 I 


some cases two sets of data had only a few lines in common for the 
determination of the reduction factors. Hence systematic errors 
with respect to both intensity and wave length may be present. It 
may be possible to detect these errors by observations of emission 
spectra, since for such spectra a larger range of intensity can be 
measured on a single spectrogram. Furnace emission spectra will 
also extend the observations to lines of higher excitation potential. 
The absorption spectrum method, while more favorable for the 
measurement of intensities of lines arising from the lowest levels, is 
limited to such lines. 


Our thanks are due to Miss Dorothy Carlson, who made a large 
number of the microphotometer tracings. 
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TEMPERATURE IN THE SOLAR REVERSING LAYER 
DERIVED FROM TITANIUM LINES 


ROBERT B. KING 


ABSTRACT 


A temperature for the solar reversing layer is derived by a comparison of the equiva- 
lent widths of solar 771 lines with laboratory gf-values. Boltzmann factors for lines 
arising from different low-energy levels are found by plotting solar equivalent widths 


as a function of gf. These give a mean temperature of 4400° K. 
An empirical curve of growth for the sun is constructed with the aid of the tempera- 


ture thus found and the laboratory gf-values for Fe 1 and 77 1. 


Comparison of the data on relative gf-values for lines of Fe1 
and 7i1 given in the preceding paper’ with the equivalent widths of 
these lines in the solar spectrum leads to a temperature for the solar 
reversing layer. Extensive measures of equivalent widths of solar 
absorption lines in the region AA 4036-6600 have been made by 
C. W. Allen,? including many of the lines whose relative gf-values 
have been determined in the laboratory. 

The problem of determining a temperature for the reversing layer 
from these data is simply a matter of finding the Boltzmann factor 
to be applied to the laboratory data which will bring the relative 
intensities of the lines arising from different energy levels into agree- 
ment with their relative intensities observed in the solar spectrum. 
| The relative numbers of atoms in two levels of energies, E, and E,, at 
a temperature T is given by the Boltzmann formula 


| —E,/KT 
(1) 


N, g,? 


| where g, and g, are the statistical weights of the levels. Taking only 
the temperature factor into account, we have 


N, (2) 


'R. B. King and A. S. King, Mt. W. Contr., No. 581; Ap. J., 87, 24, 1938. 
2 Mem. Commonwealth Solar Obs., Canberra, Australia, No. 5, 1935. 
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where Av is the wave-number difference between the two levels. 
This can be written in the form (logarithm to base ro) 


log R= (3) 


where RX is the Boltzmann factor and log R is to be added to values 
of log gf to allow for a temperature distribution of atoms among the 
various low levels from which the absorption lines arise. If solar 
values of log EW /)? are plotted against relative values of log gf for 
these lines, then all the lines arising from the same low level should 
fall on the same curve of growth. The displacement of the curves of 
growth for different low levels will be a measure of log R, and hence 
of the temperature in the part of the atmosphere where the lines are 
formed. It must be assumed, of course, that conditions there ap- 
proximate thermal equilibrium, if the temperatures derived from 
Boltzmann factors are to have meaning. 

Although laboratory and solar data for lines of both Fer and 
Ti 1 are available, only those for 771 are suitable for this purpose. 
The Fe lines, unfortunately, lie for the most part on the transition 
portion of the solar curve of growth, between the linear and radia- 
tion-damping portions, in a region dominated by Doppler effect. 
Since here EW « WV Nf + C, the displacements of curves of growth 
for lines arising from different levels along the log gf co-ordinate can- 
not be determined with any degree of precision. The Ti data, how- 
ever, are well suited for this purpose. The lines are all faint in the 
sun, and many lie on the linear portion of the curves of growth, 
where displacements are readily measured. 

Figure 1 shows a plot of log EW /) for the sun against laboratory 
values of log gf for lines of 771. The values of log gf have been re- 
duced to the lowest sublevel of the low state from which they arise 
for a temperature of 4500° K. That is, all the lines arising from the 
b’F, ;,, sublevels, for example, have been reduced to the b’F, sub- 
level, which is 11,531 wave-numbers above the o.o level a3F,. All 
lines from the b’F levels should therefore fall on the same curve of 
growth. Lines known to be blends in the solar spectrum were ex- 


cluded. 
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There is some question? as to how the solar equivalent widths 
vary with wave length. The present data do not appear to be com- 
plete or precise enough to settle this point. Plotting equivalent 
widths as various other functions of the wave length alters the 
scatter of the points but does not reveal systematic trends clearly 
enough to warrant any conclusions. The curves of growth defined by 
such plots are very nearly the same as those in Figure 1, and the 
mean temperature derived from them is not appreciably different. 


0.0 


log 


0.0 1.0 2.0 3.0 40 


log gf 
Fic. 1.—log EW/ (Sun) as a function of log gf (Lab.) 


In Figure 1 the displacements of the curves of growth from that for 
the a°F level, measured along the horizontal dotted line, give values 
of log R. Only the curves for lines from the a’P, b°F, and a’G levels 
were regarded as suitable for the purpose. The results are shown in 


Level Av log R 
8436 —1.18 4470° K 
11531 —31.67 4310 


the table, where Av is the wave-number difference between the 
lowest sublevel and a*F,. The absolute temperature T is derived by 
means of equation (3). The resulting mean temperature is 4400° K 
with an estimated uncertainty of the order of + 100°. The dotted 
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line on the right in Figure 1 gives an indication of the degree of pre- 


-cision. This would be the position of the curve for a’G relative to 


that for a’F if the temperature were 4000° K. Since the laboratory 
gf-values depend on measures with an optical pyrometer of furnace 
temperatures of 1800°-2600° C, the solar temperature derived here 
depends largely on these measures. Any error in laboratory tempera- 
ture scale and any possible departure from thermal equilibrium in 
the furnace enter into the results. 


log EW /»? 


-3.0) 
0.0 1.0 2.0 3.0 40 5.0 


log Ngf 


Fic. 2.—Curve of growth for the Sun 


It is noteworthy that no assumption need be made as to the form 
of the solar curve of growth, since the temperatures were derived 
from the displacements of empirical curves. 

The laboratory and solar data, together with the temperature 
derived here, enable one to construct an empirical curve of growth 
for the sun. This curve is shown in Figure 2, where values of log 
EW/\ (reduced to \ 5000) are plotted against laboratory values of 
log Ngf, a temperature of 4400° K being used to determine relative 
values of NV. The data for Fe and Ti were plotted separately and 
superposed as accurately as possible by shifting the two plots hori- 
zontally. The full-line curve is a theoretical curve of growth for the 
sun used by Allen. It is unfortunate that neither the Fe nor the 77 
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data extend over the whole course of the curve from linear to damp- 
ing portions. Hence a good check on the value of the damping con- 
stant cannot be obtained. There is no indication, however, of wide 
departure (except on the linear portion) from Allen’s theoretical 
curve calculated on the assumption of a value ten times the classical 
value of this constant. 

Measurements of the temperature of the solar reversing layer 
by means of spectroscopic relations involving the temperature have 
been few because of the lack of data for both laboratory and solar- 
line intensities. Real differences in such determinations might be 
expected if the spectra of vapors which absorb most effectively at 
very different atmospheric levels were examined. The temperature 
derived by Birge (4300° + 500° K)3 and that of Richardson (5700° 
+1600° K)4 from the intensity distributions in the solar bands of 
CN d 3883 and C, dA 5165, respectively, may be mentioned in con- 
nection with the foregoing result derived from neutral titanium lines. 

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 


CAMBRIDGE, MASSACHUSETTS 
September 1937 


3 Ap. J., 55,273, 1922. 
4 Mt. W. Contr., No. 422; Ap. J., '73, 216, 1931. 
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THE EFFECT OF AN ADIABATIC LAYER 
UPON SOLAR LIMB DARKENING 


PHILIP C. KEENAN 


ABSTRACT 


On the assumption that there exists a zone of instability characterized by an adia- 
batic temperature gradient beginning at optical depth 7: below the surface of the photo- 
sphere, the boundary conditions determine the distribution of temperature within the 
unstable layer. The limb darkening can then be determined for different values of 7: 
and computations carried out for boundaries of 1.0, 0.67, 0.50, and 0.28 predict curves 
of darkening similar in character to the observed curve. However, an appreciable 
effect due to the adiabatic layer could occur only if its upper boundary were considera- 
bly higher than the optical depth of two now generally assigned to it. 

According to the criterion of stability derived originally by 
Schwarzschild,’ an atmosphere in radiative equilibrium will be me- 


chanically stable if its temperature gradient, 


(ii) 
dh rad ‘ 


is less than that in an adiabatic atmosphere. Unséld? has shown 
that on the assumption of an opacity law of the form x ~ p”, the 
condition for stability corresponds to the inequality 


T+ ~ (1) 


When the value of y, the ratio of specific heats, is so low that the 
inequality is reversed, the atmosphere is unstable and convection 
currents will develop and presumably transport sufficient quantities 
of heat to establish adiabatic equilibrium within the zone of insta- 
bility. Accepting for the present Siedentopf’s conclusion that the 
convection currents will transport so little energy above the zone 
of instability that radiative equilibrium will be re-established with- 
in a very short distance,’ the inequality (1) when written as an 


! Gottingen Nachrichten, p. 41, 1906. 
f. AP., 138; 1930: 
3 See the discussion by H. H. Plaskett, M.N., 96, 402, 1936. 
4A.N., 255, 157, 1935. 
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equation defines the temperature 7, at the junction of the adiabatic 
and radiative layers. 
Since 


Tt = To(1 + 371) (2) 
we have for the optical depth corresponding to the upper boundary 
of the convective layer 


In _ ay — 9) 


or 


(4) 


Tr 


defining the optical depth at which instability will be just setting 
in for a given value of y (we shall assume that y is constant in the 
adiabatic zone). The effective ratio of specific heats at any depth 
depends upon the rate at which the various gases are being ionized. 
Computations carried out by Unsold? and by Siedentopf* show that 
the zone of instability in the sun begins at t © 2 and extends inward 
to an optical depth of about 300, i.e., over a distance of several 
hundred kilometers. 

The possibility that this region of convection influences the limb 
darkening at the surface has led Plaskett’ to attempt to derive the 
solar temperature distribution directly from the accurate observa- 
tions of darkening in different parts of the spectrum made by 
Moll, Burger, and van der Bilt on the Gornergrat in Switzerland.‘ 
Taking their data for \ 5500, Plaskett employed an approximate 
summation to solve the integral equation for B,(7) in terms of the 
intensity of the emergent radiation. He found a distortion of the 
radiative temperature distribution in the form of a nearly isothermal 
region near t = 1, followed by a gradient steeper than the radiative 
near the surface (see Fig. 7 of his paper). From this result Plaskett 
concluded that the observations indicated the presence of an un- 
stable layer of the sort predicted. 

However, the deduced temperature distribution is inconsistent 
with that which must prevail in such an adiabatic layer, where the 


5 B.A.N., 3, 83, 1925. 
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temperature gradient is at all points less than the radiative gradient. 
Furthermore, although this is undoubtedly an accurate solution 
of the equations, it is by no means obvious that it is the only tempera- 
ture distribution which will fit the data within the quite appreciable 
observational uncertainty. 

Some of this ambiguity may be removed by approaching the 
problem from a different direction, in accordance with a suggestion 
made by Chandrasekhar. Consider an atmosphere containing a 
layer in adiabatic equilibrium with its upper boundary at a specified 
optical depth 7, and extending inward to a depth so great that the 
matter beyond it will not contribute appreciably to the surface 
radiation. What will be the observed limb darkening for different 
values of 7,? 

There will be one definite curve of darkening for each 7, if the 
temperature distribution within the convective layer is uniquely 
determined by 7,, and we shall see that this is actually the case. 

Starting with the equation of hydrostatic equilibrium, 


dp = gp dx (x increasing inward) , (5) 


and the defining equation for r, 


dr = xp dx, (6) 
we have 
dr =" dp. 
(7) 
Introducing an opacity law of the form 
= Kop" (8) 
and integrating, we obtain 
n+1 
or 
n+1 
p= (r+ (9) 


where D is a constant of integration. 
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We now introduce the adiabatic condition in the form 


y-! 
T=Cip”. (10) 


Substituting from (9), and writing By(7) for the emissivity within 
the adiabatic region, we find 


4(y—-1) 


(n +1) 


Ko 


= CBx(0)[7 + (11) 


where g = 4(y — 1)/(1 + 1)y and B(o), the emissivity at 7 = 0, is 
introduced into the multiplying constant for convenience in com- 


putation. 
For a specified 7, the constants C and D are determined by the 


following boundary conditions: 
1. At 7,, By, (7;) = Bg (r:), 


where 
Br(r) = Br(o)[1 + 37]. (12) 


Combining equations (11) and (12), the condition becomes 


C= 


2. At 7, the net flux must be equal to the observed net flux (/’) 
at the surface of the sun. 


F_ 


The inward flux due to the upper layer in radiative equilibrium 
can be readily determined from the equation® 


$7:Br(0) — 2Br(o) cos sin — 3 cos p)dy . 


6 See Milne, Handbuch d. Ap., 3, 128, eq. (151), 1930. 


| | 
| 


SOLAR LIMB DARKENING 49 


Using the substitution y = 7, sec y, the solution comes out in terms 
of the exponential integral as 


3 


In order to evaluate the outward stream of radiation at 7,, we 
assume that y is constant within the convective layer. Then 


co 
0) = e(™—7) sec Ba(r) sec 0+ dr. (16) 


Substituting for B,(r7,) from (11) and writing r = t cos @ — D, 
T4(11, 0) = CBr(o)e™* %(cos + 1) — + 1)] 
and 


= 2CBx(o) { (cos see + 1) — + 1)]d(cos 0), (17) 


where I'p (g + 1) is the incomplete gamma function’ 


(71 +D) sec 0 
{ . 
° 


Solution of (17) by graphical integration is convenient and accu- 
rate, since the integrand decreases almost linearly with cos 6. If we 
write sec 6 = x and 7, + D = a, it may be shown that*® 


= 


lim 

Since the unknown D is involved in the integral, it is necessary 

to determine the constants by successive approximations. However, 

three integrations have usually been found to give sufficiently accu- 
rate values. 

With the constants known the complete temperature distribution 

is fixed, and the limb darkening can be found directly. Adding the 


7K. Pearson, Tables of the Incomplete '-Function, Cambridge, 1934. 
8 See Whittaker and Watson, Modern Analysis (3d ed., 1920), p. 341. 


50° PHILIP C. KEENAN 


contributions from the two layers, we have for the emergent radia- 
tion 


I(o, 0) = Br(r)e—* 9 sec 0+ dr + sec 0+ dr. 


Substituting for Br(r) and B,(r) and integrating, 


I(o, 0) = Br(o)[r + 3 cos — e~™ 8 + 37, + 3 cos 
+ Br(o) e? + 1) — To(q4+1)]. (18) 


For 7, = » the equation reduces to the cosine law of radiative 
equilibrium, while for 7, = o we obtain the well-known result that 
an adiabatic atmosphere shows complete darkening at the limb. 

It is important to note that the darkening for a given value of 7, 
is independent of the coefficient in the opacity law. This follows 
from the fact that » and \ enter the expressions for the constants 
only through gq, and by equation (3) 


- (10) 


for all values of . 

By means of the foregoing equations the solar limb darkening 
that would result from an adiabatic layer beginning at any specified 
optical depth can be determined. Computations have been carried 
out for 7, = 1.0,0.67,ando0.28inturn. Asan example of the numer- 
ical values, for r; = 1.0; g =0.6, D = —o.18, and C = 2.82. 

The resulting curves of darkening for the three cases are plotted 
as the broken lines in Figure 1. The straight line shows the normal 
cosine darkening for an atmosphere in complete radiative equilib- 
rium. For comparison the observations of Abbott, Fowle, and 
Aldrich as reduced to integrated light by Lindblad? are represented 
on the diagram by filled circles. 

Since the diminished temperature gradient in the adiabatic layer 
lying below the surface weakens the intensity of the light observed 
at the center of the disk but does not influence the radiation from 


9 Nova Acta Reg. Soc. Sc. Upsal., 6, No. 1, 17, 1923. 
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the limb, the effect of the unstable region is to reduce limb darkening 
over the greater part of the sun’s radius and to introduce curvature 
of the sort actually observed. However, it is clear from the diagram 
that in order to influence appreciably the observations the zone of 
instability must extend upward to optical depths less than 1. This 
last point might have been anticipated, since the contribution from 
lower depths is slight, and there would be no reason to expect any 
violent discontinuity in the temperature gradient where the two 
layers merge. 


I(0,0) 


| | | | | l | | | 
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 
cos 0 


Fic. 1 


On the other hand, we have seen that the predicted unstable 
layer is confined to greater depths (below 7 2). For this reason 
it is not possible at the present time to conclude from observations 
of limb darkening whether the prediction is correct. If later there 
develop grounds to believe that the region of instability can affect 
the surface layers, it will be worth while to refine our assumptions 
in order to predict exact temperature gradients. The general char- 
acter of the curves of darkening to be expected has been established 
by our present investigation and is not altered greatly by even con- 
siderable changes in the constants. 
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It is of incidental interest to see how closely the observations 
could be fitted if we were permitted to choose an appropriate value 
of 7,. Figure 2 shows a curve of darkening for 7, = 0.5, for which a 
weighted mean value of y within the adiabatic layer was used to 
improve the approximation. The computation was carried out for 


I(0,0) 15500 

e 
Z 
0.8 
a 

= 
0.4 

l | l | | | | 


0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 
cos 6 


FIc. 2 


5500 by an approximate summation, using Plaskett’s convenient 
tables, in order to allow comparison with the measures in this 
spectral region by Moll, Burger, and van der Bilt, indicated by the 
filled circles in the diagram. The fit is close except for cos 8 = 0.4, 
where the observations may be less certain because of atmospheric 
scattering and scintillation. 


In conclusion I should like to acknowledge many helpful dis- 
cussions of this problem with Dr. S. Chandrasekhar. 


YERKES OBSERVATORY 
September 1937 
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LINE INTENSITIES IN ULTRAVIOLET 
STELLAR SPECTRA 


LAWRENCE H. ALLER AND LAURENCE G. STODDARD 


ABSTRACT 


Eye-estimates of intensity, on an arbitrary scale, of 17 ultraviolet lines of neutral 
iron, magnesium, and aluminum have been made on 606 plates distributed among 31 
representative supergiants, giants, and dwarfs. These intensities are compared with 
half-breadths computed on the basis of Pannekoek’s recent theory, and the results are 


discussed. 

In a recent paper Pannekoek' has attempted the theoretical calcu- 
lation of line intensities in stellar spectra. Starting with an assumed 
composition of a stellar atmosphere based on Russell’s work,? and 
using formulae for the stellar absorption coefficient supplied by Men- 
zel and Pekeris,* he has deduced numerical values for the coefficients 
of opacity and absorption. These data have been utilized in his 
numerical integration of the equations of radiative transfer, making 
possible the theoretical determination of the half-breadths of spectral 
lines. 

The present writers have undertaken an observational test of this 
theory, using eye-estimates of intensity on an arbitrary scale. For 
this purpose a series of plates, taken by Professor C. D. Shane with 
the quartz spectrograph attached to the Crossley reflector of the 
Lick Observatory, was used. While the theory gives half-breadths 
and observations yield eye-estimates of intensity, certain features 
can be examined, such as the theoretical predictions of maxima or 
the complete absence of lines. 


THE OBSERVATIONAL MATERIAL 


The observational material consists of 172 spectrograms taken 
during the summer of 1926. These plates were distributed among 
65 stars. Of this group, 27 stars were rejected for one of the following 
reasons: (1) nebulous lines or definite evidence of rotation; (2) 
spectroscopic binaries with two spectra visible; (3) poor plates. 

* Pub. Astron. Inst. Amsterdam, 4, 1935, and Addendum. 

2 Ap. J., 70, 11, 1929; 75, 337, 1932; and 78, 281, 1933. 3M.N., 96, 77, 1936. 

53 


54 LAWRENCE H. ALLER AND LAURENCE G. STODDARD 


Under (1) a number of early-type stars were rejected which showed 
fuzzy lines or possessed a high rotational velocity (greater than 100 
km/sec), according to the studies of Struve and Elvey and Miss 


TABLE 1 
Spectral No. of 
Star Vis. m Vis. M 

e Serpentis........ A6s 3.8 1.8 2 
y Virginis.......... Fo 3.8 2.8 3 
cF4 —0.9 4 
Fs 4.0 4.0 3 
¥ Serpentis......... F5 3-9 2 
F6 4.1 2:6 4 
ec F7 2.8 3.2 2 

d Serpentis......... Gr 4.5 Act 2 
B Draconis......... cG2 3.0 —1.7 2 
37 Leo Minoris...... cG2 4.8 —1.6 I 
G2 3.4 0.6 2 
cG3 —1.3 I 
G4 3-8 0.8 I 
G4 3.6 0.3 2 
pe G4 3.4 4.2 2 
a Ursa Majoris..... G6 2.0 0.0 3 
o Draconis......... Go 4.8 5.6 I 
Ko 2.7 —o.I I 
0.2 0.2 4 
x Ursa Majoris..... Kr 3.8 0.4 2 
cKr 4.0 —1.3 2 
70 Ophiuchi......... Kr 4.3 5.9 2 
8 Ursa Minoris..... K5 2.2 —0.5 I 
+ Draconis......... Ks5 2.4 —o.I 3 
K6 5.6 I 
8B Andromedae...... Mo 2.4 0.2 2 
Mr 1.2 —3.8 I 
M2 2.6 —o.6 2 


Westgate.’ It was necessary to select stars without appreciable rota- 
tional velocities, since this theory does not take account of rotation. 
Since the lists of Struve and Elvey and Miss Westgate are not com- 
plete, it was possible for us to reject only those stars which their lists 


4 Ap. J., 7%, 227, 1930; 77, 141, 1933. 
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had also contained. However, the stars with nebulous lines were all 
rejected, as this phenomenon was believed, in general, to be due to 
rotation. Spectroscopic binaries with two spectra visible were re- 
jected for obvious reasons. 

For some stars only one or two plates were available, and these 
were sometimes either underexposed or did not have excellent defini- 
tion. Among the 38 remaining stars, seven were of early type and 
were not used in the present investigation, which deals with lines of 
neutral iron, magnesium, and aluminum. The hot stars usually con- 
tained too few lines of any kind. In the final study, 66 spectrograms 
of highest quality were used. The plates were distributed among a 
representative set of giants and dwarfs, with a few supergiants (see 
Table 1).5 

THE SELECTION OF THE LINES STUDIED 

The spectra cover the range from \ 3500 to the H and K lines. In 
the cooler stars the number of lines is large, and the blends are very 
numerous. The most prominent lines are due to iron with fainter 
ones of titanium, nickel, cobalt, aluminum, and magnesium. Plates 
of Arcturus, the sun, and x Draconis were measured and reduced by 
the Hartmann formula, on the basis of the iron comparison spec- 
trum. The Revised Rowland Table was used for identifying the lines. 
The only lines retained for purposes of comparison were either un- 
blended or possessed one predominant component, preferably 
blended with a line of the same element in the same state of ioniza- 
tion. Table 2 gives a summary of the lines studied with the multiplet 
designations and Rowland intensities. Some lines beyond X 3500 
were studied; but the weakness of the continuous background 
vitiated the estimates, and a discussion of these lines is not included 
here. 

ESTIMATES OF THE INTENSITIES 

The good plates were divided into six classes, according to the 
nature of the continuous spectrum. The lines to be studied were 
marked on an enlargement of the spectrum of the moon. The plates 
of each group were thoroughly mixed, and the intensities of the 
chosen lines were estimated visually on a scale of 15. Estimates of 
the intensities of all the groups were carried out, and the procedure 


5 Ibid., 81, 187, 1935; Mt. W. Contr., No. 511. 
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was repeated until the intensities on each plate had been estimated 
five times. Means were then computed. 

The reason for dividing the plates into these groups and estimat- 
ing the intensities for each group separately was to see if there was 


TABLE 2* 
Ground 
Line Energy Rowland J Multiplet 
(Volts) 
Aluminum: 
A 3044.02..... 0.000 15 3?P°— 
3901.54..... 0.014 20 3°P°— 47S 
Magnesium: 
A 3838.50..... 2.705 2-1 25 33P°— 33D 
Tron: 

0.855 5-4 25 a5F —ysD° 
ISO 0.gII 4-3 20 a5F — 
3849.98..... 1.007 I-o 10 a5F — ysD° 
3734.88..... 0.855 5-5 40 — ysF° 
3740. 50...+. 0.gII 4-4 20 — ysF° 
9763.80..... 0.986 2-2 10 a5F — ysF° 
| 1.007 I-I 8 aSF — ysF° 
3709.26..... 4-3 8 — y5F° 
0.000 4-5 40 asD —25F° 
0.051 3-4 30 —z5F° 
0.855 5-6 30 —25G° 
3008.87... 1.007 I-2 20 asF —25G° 
2037. 0 0.954 3-4 15 asF —z5G° 
SEBO, 0.gII 4-5 20 —73G° 


* From Miss Moore’s A Multiplet Table of Astrophysical Interest. 


any dependence of the intensity estimates on the density of the 
background. It was found that, when the results were plotted, there 
was not much difference, except at the ends of the spectral range, for 
plates of the extreme groups. The plates were all Eastman 4o’s, 
except for one Eastman 33. 
THE THEORETICAL CALCULATION OF THE INTENSITIES 
ON THE BASIS OF THE PANNEKOEK THEORY 

If we confine our attention to lines sufficiently strong, so that the 
point at which the intensity falls to half that of the continuous back- 
ground is within the resonance wings and not within the Doppler 
core, we may use Pannekoek’s abbreviated formula 


(any ( 6 Axyf 
4000/ y s ’ 
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where 6/y¥ is the ratio of the quantum mechanical to the classical 
half-breadth of the scattering coefficient, f is the oscillator strength 
related to the Einstein probability coefficient, A is the abundance, x 
is the percentage in the proper state of ionization, y is the percentage 
excited to the ground state of the transition in question, and s is the 
appropriate scattering coefficient. In practice we compare the same 
line throughout the spectral range. Thus \/ 4000, 6/y, and f are con- 
stant and A is assumed constant, so that for a given line we may put 
(for two stars of temperatures, 7, and 7,) 


(An), 


The percentage of atoms in a given state of excitation and ioniza- 
tion was computed by means of the Saha and Boltzmann expres- 
sions. The partition functions were computed with the aid of the 
data of Bacher and Goudsmit. Modifications of the Boltzmann 
formula for deviations from thermodynamic equilibrium were not 
considered. The effect of the dilution of the radiation upon the 
ionization was taken into account. 

Pannekoek considers the problem of line formation under the fol- 
lowing three conditions: 

Case A.—The number of atoms is decreasing, owing to ionization, 
as the temperature increases (the case of neutral atoms in the hotter 
stars). 

Case B.—The number of atoms is constant or a maximum (e.g., 
the case of neutral atoms of metals such as Fe, Cr, Nz, etc., at low 
temperatures, or of Ca before the second ionization sets in but when 
the neutral calcium is all ionized). 

Case C.—The number of atoms is increasing, owing to ionization. 
This case is of significance only for the ionized atoms (see Fig. 1). 

Now s is the value of the scattering coefficient for a point in the 
line for which the intensity falls to one-half that of the continuous 
spectrum. When all the atoms are in the proper state of ionization, 
s is equal to s,. We may write for Case A, 


¥ 
P+ K’ 


5S = So 
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by virtue of the well-known ionization formula (aP/|1 — «]) = K. 
The electron pressure P is much smaller than K for the case of high 
ionization, so that for the purposes of his numerical integrations 
Pannekoek took s = (s.P)/K = oP, where o = s,/K. For Case B 
he merely took s as constant. In the present investigation we need 
concern ourselves only with Cases A and B, since we are dealing 
with lines of neutral atoms which remain approximately constant in 


Fic. 1.—Pannekoek’s three cases for a hypothetical element X and X*, e.g., Mg and 


number over a certain temperature range and then decrease for 
higher temperatures. 

Pannekoek tabulates o and s with the surface gravity g, 5040/7, 
and the wave length as arguments. Now in Case A, s, was taken as 
equal to oK in the numerical integrations. Actually, we should have 
had (K + P) instead of K, so that, after having computed K and y 
and after having taken o and P from the tables, we have 


For the Case B we look up s in Pannekoek’s Table 29a, and, noting 
that s = (s.P)/(P + K), 


(Ad)? 


For the purposes of the numerical integrations, the neglect of these 
proper factors was justifiable; but in comparing a line from star to 
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star where the surface gravity and temperature vary (and hence K 
and P must vary), we must include them, to get the best results. 
Usually, the intensity of a line was computed on the basis of both 
cases, A and B, for those intervening regions in which neither case 
is strictly valid. 

We need to know the surface gravity and temperature for each 
spectral class among the stars studied. The surface gravity and 
temperature are connected by the empirical formulae given by 
Russell :° 


logio =—0.65+0.285 (giant sequence), 


logis = + 0.76 — 0.972 (main sequence). 


For the sun, log g, is 4.4; 3 = 11,600/T. 

The temperature scale adopted was that given in Russell, Dugan, 
and Stewart’s Astronomy.’ For use in the ionization formula 
1 —« = P/(P + K), the mean electron pressures, P,,, are taken 
from Pannekoek’s Table 22, and K is computed in the usual manner. 
The o and s are found by interpolation in Tables 28 and 29. For the 
supergiant stars no temperature scale was available, and we have 
assumed the following values by the extrapolation of the difference 
between giants and dwarfs for the temperature of a given spectral 
class: 

Fo Fs Go Gs Ko Kg5 M2 Ms 
7400 6500 5500 4200 3700 3200 3000 2600 

On this basis the absolute bolometric magnitudes for the stars 

used were computed from the expression 


M, — M, = 10 log T + 91500 — 42.1, 


given by Russell, Dugan, and Stewart. Use was then made of a 
table published by Mrs. Payne-Gaposchkin in her Stars of High 


6 Ap. J., 78, 270, 1933; Mt. W. Contr., No. 477, 1932. 
TP. 
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Luminosity,’ to find appropriate surface gravities for the supergiant 
stars. We adopt the following: 


Temperature.... 7200 6300 5600 5040 4200 3600 3150 2520 
3.0 2.7 2.5 2.2 1.8 1.0 


COMPARISON OF THEORY AND OBSERVATION 


The eye-estimates of intensity were plotted against temperature 
for each line in the dwarf, giant, and supergiant sequences. The best 
line through the points was drawn in each case. Sometimes the 
scatter was bad, but for the most part the observations permitted 
fairly definitive curves to be drawn. Typical curves are reproduced 
in Figure 2. 

The logarithms of the half-breadths, as computed by the method 
outlined above, were plotted against the temperature. For purposes 
of further comparison between theory and observation, the actual 
half-breadths were plotted also. A comparison of our scale with the 
Rowland scale was also made by use of the plates for the sun. In re- 
gard to the scatter of the individual values of the intensity estimates, 
it should be remarked that not all of the scatter is due to errors of 
observation. Some of it may certainly be attributed to peculiarities 
of the individual stars, which, for a given spectral class, often showed 
considerable range in absolute magnitude. 

Aluminum.—The only aluminum lines available were \ 3944 and 
d 3961. Theory predicts that above 5000” the lines should be slightly 
stronger in dwarfs and giants than in supergiants, but that below 
5000’ the lines should be strongest for the supergiants with a maxi- 
mum half-breadth of about three times that of the dwarfs. Flat 
maxima are predicted at 4000° for the dwarfs and near 3000” for the 
giants and supergiants. 

Actually the lines appear at a higher temperature in the dwarfs 
and reach a greater maximum intensity than in either the giants 
or the c stars. The maximum is poorly defined in all cases. The lines 
are strongest for the dwarfs, next strongest for the giants, and weak- 
est for the supergiants. The mean maximum intensity is 4.5 for the 
c stars and 8 for the dwarfs, on our scale. 
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Fic. 2.—Observed intensity of \ 3631, Fe and Cr*; (a) dwarf, (6) giant, (c) super 
giant. The temperature is in degrees absolute. 


| 
| 


62 LAWRENCE H. ALLER AND LAURENCE G. STODDARD 


Magnesium.—For magnesium, \ 3838, theory predicts lines of 
about equal intensity in giants, supergiants, and dwarfs, down to 
about 5000°. Below sooo” the line should be strongest in giants and 
supergiants and weakest in the dwarfs. (At 4200° the line should be 
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Fic. 3.—The intensities of the aluminum lines \ 3944 and d 3962 in dwarfs, giants, 


and supergiants; (a) observed and (6) theoretical curves. 


twice as strong in the supergiants as in the dwarfs.) The predicted 
maxima for the supergiants, giants, and dwarfs are, roughly, at 
4000°, 4100°, and 5000’, respectively. 

The observational data are rather uncertain, but they definitely 
show that the line appears at a higher temperature in the dwarfs, 
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and for the higher temperatures is much stronger in these stars than 
in the giants and supergiants. The maximum intensities appear to be 
about equal for dwarfs and supergiants; but the maximum for the 
giants seems slightly weaker, probably the effect of observational 
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Fic. 4.—The intensities of the magnesium line \ 3838; (a) observed and (6) theo- 
retical. 


errors. The maxima for giants and supergiants appear at 4000° and 
3500°, respectively; and the maximum for the dwarfs is not indi- 
cated. 

Tron.—For iron, the material is the most extensive and we may 
draw the soundest conclusions. For all the lines except those beyond 
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the Balmer limit, the theory predicts roughly equal intensities for 
giants, supergiants, and dwarfs for any temperature above 5000”. 
For the aS‘F — z'G multiplet the lines should be strongest for the 
dwarfs at the higher temperatures, but in all three cases the intensity 


Int. 


= 


gooo 7000 5000 3000 af 
Fic. 5a.—Observed intensities of the iron lines \ 3720 and \ 3737 in dwarfs, giants, 


and supergiants. 
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Fic. 5b.—Observed intensities of the iron line \ 3826 in dwarfs, giants, and super- 
giants. 


should fall rapidly as the temperature increases. The theory predicts 
strong maxima for the supergiants with half-breadths ranging from 
three to five times those of the same line in the dwarfs. For the lines for 
which the excitation potential is of the order of 1 volt, the tempera- 
tures of the maxima are, respectively, 3600°, 3900°, and 4700°, for the 
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supergiants, giants, and dwarfs. For the lines whose ground level is 
about o volts, the maxima are at 3500°, 3600°, and 4200°. 

In all cases the dwarf lines are much stronger than the giant and 
supergiant lines for the higher temperatures, and the giant lines are 
usually more intense than the supergiant lines. The maximum in- 
tensities of the lines are certainly just as great for the dwarfs as for 
the supergiants, and probably greater. The a‘D — z5F multiplet 
yields a very uncertain maximum; but for the other multiplets of 
iron studied, the positions of the maxima seem to be nearly the same 
as that predicted by theory for the supergiants. The agreement is 
generally good for the giants; the single discrepancy is probably due 
to errors of observation. For the dwarfs the maxima are near the end 
of the range of temperature to which the observations extend, and 
therefore the agreement between theory and observation cannot be 
tested. 

A bad discrepancy between theory and observation exists at the 
head of the Balmer series. Theory predicts that the lines of the 
a5F — z5G° and a5F — z3G° multiplets, AX 3581, 3570, 3609, and 
3631, should disappear near 7400°; but they persist to 8400° along 
the main sequence. Data for the giants are not available for the 
higher temperatures, but the lines are still strong for the hottest star 
in the list. This contradiction between theory and observation also 
holds for the supergiants. 

To summarize our results, we may say that the theory tends to 
predict too intense lines for the giants and supergiants, as compared 
with the dwarfs, especially at the higher temperatures. Excluding 
the case of magnesium, for which the data are rather weak, the 
maxima, as predicted by theory, agree quite well with the observa- 
tions. A notable discrepancy between theory and observation seems 
to occur for the lines below the Balmer limit. These lines should fall 
off in intensity very rapidly where hydrogen begins to be appreciably 
ionized, but they persist to much higher temperatures than one 
would expect on the basis of the theory. The discrepancy seems to be 
real, and it is not possible to account for it on the basis of blends. 


9 At higher temperatures, \ 3631.5 Fe is blended with \ 3631.7 Crt, atP — z4P°. 
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We are indebted to Professor Shane, who suggested the problem 
and furnished the observationai material and under whose guidance 
the investigation was carried out. We are also thankful to Mr. 
James Fidiam, of San Francisco, for his beautiful enlargements of 
many of our plates and to Mr. Gordon Wares, of the Yerkes Observ- 
atory, who helped in classifying the plates. 
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THE FUNDAMENTAL ARC SERIES IN SODIUM 
AND ITS PRESENCE IN THE SUN 


PAUL ROOD AND R. A. SAWYER 


ABSTRACT 

Using a powerful discharge tube of new design, the infrared arc spectrum of sodium 
has been photographed, and measured with a thermopile, as far as \ 12679, 3?7D—5?F°. 
In addition to this line, five new members of the ?D?F° series in Na I were measured 
with an accuracy of +0.2 or +0.3 A. Data supplied by Miss C. E. Moore and Mr. 
H. D. Babcock, of the Mount Wilson Observatory, indicate that this sodium series is 
present in the solar spectrum. New values for the ?7F° terms are given. In addition, 
three members of the 4?S— 5?P° series were observed for the first time. 

In connection with their program of identification of the infrared 
solar lines at the Mount Wilson Observatory, Miss Charlotte E. 
Moore and Mr. H. D. Babcock have called to our attention the fact 
that no recent or accurate measures are available of the Bergmann 
series in the arc spectrum of sodium. The approximate positions of 
the lines can be predicted within a few angstroms, but no values good 
enough for certain identification are known. 

The first measurement in the sodium ?D?F° series was made in 
1907 by Bergmann," who in his original work observed, but did not 
classify, the second member of the series. In the next year, both the 
first and second members were measured by Paschen,? who gave 
their wave lengths as \ 18459.42 + 2. and \12677.54 + 8. No 
other measurements in the fundamental series appear to have been 
made, although a calculated 6?F° term has been obtained from a 
supposed 3?P° — 6?F° combination at \ 4665 and \ 4660. This term 
does not agree well with the one determined in the present work. 

The successful measurement of the lines of the fundamental series 
requires a vacuum source, since in open arcs they appear broad and 
fuzzy, if they appear at all. The source finally found most suitable 
was the positive column in a large end-on discharge tube. The tube 
was of pyrex glass, 3 cm. in diameter and about a meter long, with 
large cylindrical aluminum electrodes in side arms. Pure redistilled 
sodium was introduced into the tube in glass boats, and the re- 


* Zs. f. Wiss. Phot., 6, 113 and 145, 1908. 2 Ann. d. Phys., 27, 543, 1908. 
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quired pressure of sodium vapor was produced by an electric furnace 
surrounding the tube. The discharge was carried by pure helium or 
neon. The series was more intense in helium; but the third and fifth 
members fell too close to helium lines for measurement, and a dis- 
charge in helium-free neon was used to obtain these lines. The heli- 
um or neon was circulated through charcoal traps in liquid air and 
through the tube by means of a glass mercury pump. Sixty-cycle 
current from a 2000-volt transformer at currents of 1 to 1.5 amperes, 
at suitable pressures of sodium and of the gas used, provided a very 
intense sodium light. 

The intensity of the sodium spectrum in the discharge was 
sufficient to permit thermopile measurements in the infrared with a 
grating spectrograph. This instrument,’ which was placed at our 


TABLE 1 
MEASUREMENT OF THE 3?P°—4?S DOUBLET OF SODIUM 
Classificati Wave Length Av Av E 
Observed Observed Calculated 


disposal by Professor E. F. Barker, who assisted in its use, has 
mirrors of 100 cm focal length and a r1oo lines per centimeter grat- 
ing. Using the helium line \ 10830.3 as a standard, the second mem- 
ber of the sodium fundamental series, 37D — 5’F°, was measured as 
 12679.0 + 0.5 A. The contour of this line was mapped, and its 
width was found to be not greater than the slit width of the spectro- 
graph, an indication of the sharpness of the lines from this discharge. 
As a test of the accuracy of this measurement, the 3?P° — 47S 
doublet, whose separation can be calculated from known terms, 
was measured, with the results shown in Table 1. 

The shorter members of the fundamental series were photo- 
graphed with a Littrow-type flint-glass spectrograph having a dis- 
persion, at 10,000 A, of about 70 A per millimeter. Eastman M, 
(), and Z plates were used in the regions of their respective optimum 
sensitivities. The exposures ran from 3} to 5 hours. The helium 


3 Barker and Meyer, Trans. Faraday Soc., 25, 912, 1929. 
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or neon lines emitted by the discharge tube, depending on the gas 
used, were taken as standards, the values being those from the 
recent measurements at the Bureau of Standards.* For converting 
the wave lengths to wave numbers, the values for the dispersion 
of air by Meggers and Peters,’ extrapolated above \ 10,000, were 
used. 

The final values adopted for the fundamental series lines of 
sodium are given in Table 2. The *F° term values and effective quan- 
tum numbers are also given. 


TABLE 2 


THE FUNDAMENTAL SERIES OF SODIUM, 3?D—m?F° 
37D = 12276.18 


Term Wave Length |Wave Number m?ke Neff 
18459. 5 5415.8 6860.4 3.9904 
6°F° 10834.4 9227.3 3048.9 5.9903 
9961.0 10030.4 2239.8 6.9905 
9466.0 10561. 2 1715.0 7.9991 
9154.7 10920.4 1355.8 8.9905 

8943.6 11178.1 1098.1 9.9966 


* Paschen’s value, not measured in the present work. 
¢ Calculated term. Interference in both helium and neon. 


With the exception of \ 8796, which appeared on only one plate, 
the lines were each measured on two to four plates. The consistency 
of the various values indicated a reliability of +o.2 or +0.3 A 
for these lines. A further direct check on the accuracy was given by 
measurements on three members of the 47S — m?P° series, whose 
wave lengths can be computed from their known terms, determined 
from very accurate measurements at shorter wave lengths. These 
lines were here observed for the first time. The comparison is given 
in Table 3. 

Miss Moore and Mr. Babcock have supplied us with information 

4 Meggers and Dieke, Nat. Bur. Stand. J. Res., 9, 121, 1932, Meggers and Humphreys, 


ibid., 10, 427, 1933. 
5 Bull. Bur. Stand., 14, 698, 1919. See also Babcock (Phys. Rev., 46, 382, 1934) ona 
method for obtaining infrared wave numbers from Kayser’s Tabelle der Schwingungs- 


zahlen. 
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on the closest solar line to each measured sodium line of the funda- 
mental series and also its probable solar identification. This infor- 
mation is given in Table 4. It is evident that this sodium series 
appears in the solar spectrum. 


TABLE 3 
NEW LINES IN 4?S—m’P° SERIES OF SODIUM 


Classification Observed Calculated* Error 
7809 . 4 7809 .6 


* The calculated value is the mean wave length of the doublet. 


TABLE 4 
THE FUNDAMENTAL SERIES OF SODIUM IN THE SUN 


LABORATORY SUN 
TERM 
Wave Nearest Int.* Solar ¥ Solar 
Length Solar A _ minus Lab. X Ident. 
18459.5 | Too far out 
12679.0 80.40 ON +1.4 Na 
10834.4 34.01 5 —0.4 Na 
9961.0 61.38 ° +0.4 Na 
9406.0 65.98 —3? —0.0 Na 
_122 —1N) Atm. 
55-4 Atm. 
8943.6 43.058 —0.5 Fe 
8796. 96.491 +0.5 Fe 


* Underline denotes that estimated intensity is from more than one plate. 
+ Paschen value of wave length. 
t Wave length calculated. 


We are much indebted to Dr. H. G. Beutler for aid in the distilla- 
tion of the sodium and in the design of the tuBe, to Professor 
E. F. Barker for his co-operation in the infrared grating measure- 
ments, and to Miss C. E. Moore and Mr. H. D. Babcock, who 
instigated this work and supplied the solar data. 


DEPARTMENT OF PuysIcs 
UNIVERSITY OF MICHIGAN 
October 1937 


| | 

| | 
| 
| 
| 
| 
t 
{ 
| 


NOTES 


THE ORBIT OF b PERSEI 


ABSTRACT 
A new orbit determination of b Persei confirms Harper’s orbit and supports his view 
that changes are taking place in some of the elements. 
Previous orbit determinations of b Persei (a[tgoo] = 4'10™7, 
[1900] = 50°03’, visual magnitude 4.57, type A2) have been made 
by Cannon’ from Ottawa observations of 1910-1911 and by Harper’ 


TABLE 1 
Date (G.C.T.) Phase Velocity | o-c (Final) 
(Km/Sec) 

1934 Sept. 26.306.... ©. 000 — 9.2 +3.8 
0.521 +58.6 —4.2 
Oct. 25408... 1.517 —10.7 +3.5 
0.3600 +55-9 +6.7 
SAL? 1.414 — 23.0 —o.8 
0.205 +32.5 —I.1 
29) 2017 0.026 =10.4 
20-401... 0.493 +63.9 +2.2 
Nov. 1.378.... 0.942 +14.8 —7.4 
0.870 +41.9 +8.3 
24.432....] 1.086 — 1.8 5 
Dec. 0.163 +12.7 —I.1 
1.159 — 8.7 +o0.8 
1935 Jan. 5.354.... 0.241 +29.2 +o.1 
0.281 +36.8 +0.3 
0.108 + 2.2 —1.4 
0.133 + 0.6 —7.6 
0.162 +-17.1 +3.3 
20:000...; 0.494 +58.4 —3.3 


from Victoria observations of 1929-1930. Harper’s orbit differed 
from Cannon’s in such a way as to indicate real changes in some of 
the elements. 

The material for the present orbit determination consisted of a 
series of 19 spectrograms made with the Bruce single-prism spectro- 

™ Pub. Dom. Obs. Ottawa, 1, 283, 1914. 

2 Pub. Dom. Ap. Obs. Victoria, 4, 309, 1930. 
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graph attached to the Yerkes 40-inch telescope, the series extending 
from September 26, 1934, to January 26, 1935. On these spectro- 
grams the lines measured were AA 3933.27, 4101.74, 4340.47, 4481.23, 
and 4549.53- 

The period adopted by Harper, 1.52738 days, fitted the observa- 
tions better than any other; accordingly, it was used. The data for 
the individual observations are shown in Table tr. 


Km/sec 


Ook ° 


20F 


—20P 


0.0 0.8 1.2 1.6 Days 


Fic. 1.—Velocity-curve of b Persei 


Four observations were put into two groups, resulting in 17 
normal places. Preliminary elements were chosen, using a graphical 
method by R. K. Young. They are as follows: 


P=1.52738 days K=43.5 km/sec 
€=0.05 +20.6 km/sec 
w= 270° T=JD 2418956.643 


Since the value of e was so small, 7 was fixed and a least-squares 
solution was carried out for corrections to e, w, K, and y, resulting 
in small corrections to the preliminary elements. The final elements, 
together with the computed probable errors, are shown in Table 2, 
along with Cannon’s and Harper’s elements for comparison. The 
probable error of a single observation is: for Cannon’s observations, 
+ 3.9; for Harper’s, + 2.4; and for the writer’s, + 2.8. 

Harper has pointed out that the changes in y and K between 
1910-1911 and 1929-1930 are probably real. The 1934-1935 orbit 
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indicates the probability of further real changes in y and K since 
1930. The change in K is 4.2 km/sec; the change in y is 6.3 km/sec; 
but, applying the corrections suggested by Moore? to the Victoria 
and Yerkes velocities, the change in y is reduced to 4.3 km/sec. 
These changes are well in excess of the probable errors and are prob- 


TABLE 2 
Cannon Harper Yerkes 
Element 
(1910-1911) (1929-1930) (1934-1935) 

1.52732 days 1.52738 days 1.52738 days 
0.22 0.017 0.047 +0.030 
I51°75 271°53 270°95+1°30 
+23.09 km/sec +13.53 km/sec +19.83+0.86 km/sec 
ee er 41.89 km/sec 38.85 km/sec 43.10+1.37 km/sec 
JD 2418956. 166 JD 8956. 662 JD 8956. 643 
837,600 km 815,840 km 805,440 km 


ably real. No real changes are indicated in the other elements since 
1930. The period of 1.52738 days connects the 1929-1930 and the 
1934-1935 observations as well as could be expected in view of the 
changes in y and K. 

Although Cannon measured a second spectrum, neither Harper 
nor the writer found any evidence of it on the Victoria and Yerkes 
plates, respectively. 


The writer wishes to express his thanks to Dr. Otto Struve for 
suggesting this orbit determination, the measures for which the 
writer made at the Yerkes Observatory in 1935 during the tenure 


of an 1851 Exhibition Scholarship. 
J. F. HEARD 
Davin DUNLAP OBSERVATORY 
September 1937 


NOTE ON THE GLOBULAR CORONA 


The two accompanying eclipse photographs were made at Moro, 
Peru (long. +78°10/5; lat. —9°8/5; elevation 2400 ft.), on June 8, 
1937. The first (a), here reduced in size, was made with a 5.5-inch 
visual objective 72 inches in focal length, belonging to Swarthmore 
College. It represents an 8-second exposure in a scheduled series 


3 Pub. Lick Obs., 18, xii, 1932. 


PLATE I 


THE SOLAR CORONA OF JUNE 8, 1937 
(a) Exposure, 8 sec., {/13; (b) exposure, 2 sec., f/3.5 
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made on Eastman 1G plates with the aid of a clock-driven polar 
axis. 

The second photograph (0), enlarged to the same scale, was made 
with a Retina camera equipped with a lens 5 cm in focal length and 
operating at {/3.5. This picture, one of two made with the Retina 
(mounted on the same polar axis), was exposed 2 seconds, the other 
4 seconds. Both the Retina pictures were made on Eastman day- 
light Kodachrome and were processed in the usual manner at 
Rochester. 

Picture (6) shows several things. (1) The extent of the halation 
around the moon’s disk is apparent when the two pictures (a) and 
(b) are compared. The disk of the moon is a dull gray, best described 
as mouse-colored. (2) Picture (b) shows the outer corona ordinarily 
photographed. On the original color-slide this is badly overexposed 
but is pale cream yellow, almost colorless. (3) This picture shows 
the globular corona inclosing the irregular outer corona. 

This is a more adequately exposed photograph of the globular 
corona noted first by Major A. W. Stevens in some of his aerial 
photographs of the same eclipse made at an elevation of 25,000 
feet off the coast of Peru. The reality of this globular envelope was 
tested on Stevens’ films at the University of Rochester’s Institute 
of Optics and in the Eastman Research Laboratories, where it was 
definitely shown not to have been introduced by the apparatus.’ On 
the color-slide this ring stands out clearly as a faint sand-colored 
disk contrasting sharply with the greenish blue of the sky im- 
mediately surrounding the sun. The diameter of the globular 
corona on the original slide is about 2°5. The 4-second exposure, too 
difficult to reproduce, shows a nearly identical halation-reduced 
lunar disk and an outer corona of the same extent, but shows the 
globular corona as a fainter disk 4° in diameter standing out less 
sharply against a lighter sky. Beneath the sun in this picture the 
globular corona is merged with the fainter extensions of the salmon- 
colored light around the horizon a few degrees away. The apparent 
color difference between the globular corona and the surrounding 
sky would seem to be of value in obtaining further observations. 

A photograph showing the globular corona was obtained by 


t Personal communication. 
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Putilin at the Siberian eclipse of June 19, 1936. Doubtless other 
photographs exist of this hitherto unnoticed globular corona from 
earlier eclipses. An investigation of old plates might add consider- 
ably to the meager information now at hand. Planned study of this 
corona is a possibility for the next useful eclipse, in 1940. 

The two color pictures and the series with the Swarthmore lens 
were made by me, as the representative of the Junior Astronomy 
Club of the American Museum of Natural History, on the Hayden 


Planetarium—Grace Eclipse Expedition. 
Dana K. BAILEY 
STEWARD OBSERVATORY 
Tucson, ARIZONA 
September 1, 1937 


PHOTOELECTRIC MEASURES OF HD 199140 


The probability of a variation in the light of HD 199140 was 
suggested in a paper read by Dr. R. M. Petrie at the fifty-eighth 
meeting of the American Astronomical Society. The period deduced 
by Dr. Petrie from variations in the spectrum was 0.20103 day. 

This star was first observed with the photoelectric photometer of 
the Washburn Observatory on September 25, 1937, and the varia- 
tion was immediately detected. Observations were made on ten 
nights, from September 25 to November 10. The comparison star 
was HD 198820. The positions are as follows: 


Star R.A. (1900) Dec. (1900) Mag. | Spec. 
20so™1 +28° 8’ 6.44 Bis 
20 47.9 +32 28 6.30 B5sk 


The light-curve is given in Figure 1; 111 observations were 
combined to give 22 normals. The total variation is 0.180 mag. The 
phase, in decimals of a period, is computed from the maximum by 
means of the elements 

Max. = JD 2428802.6689 + of20103 + E 
Max. — Min. = 041235. 


7 Ap. J.,85, 341, 1937. 
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The short interval of the observations does not permit a correction 
of the period. The normal magnitudes are given in Table 1. The 
rise from minimum to maximum is slower than the fall from maxi- 


A Mag.[— 


0.05} 
6 


0.0 0.2 0.4 0.6 0.8 1.0 1.2 Phase 


Fic. 1.—Normal magnitudes of HD 199140 


TABLE 1 
NORMAL MAGNITUDES OF HD 199140 


Phase A Mag. Residual Phase A Mag. Residual 
.078 — .002 (07 .129 — .o16 
. 209 — .002 .097 + .009 


mum to minimum, the former requiring 0.614 and the latter 0.386 of 
the entire period. This star belongs to the 6 Cephei class of variables 
but does not appear, from this short series of observations, to show 
any certain variation of amplitude. 

C. M. HUFFER 


WASHBURN OBSERVATORY 
UNIVERSITY OF WISCONSIN 
November 1937 
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NINE WHITE STARS OF LARGE PROPER MOTION 


The following nine stars with spectral types earlier than F5 have 
been found among some 160 objects with proper motions larger than 
o”3. The complete data for these objects will be given at a future 
occasion, when more stars will have been observed. The magnitudes 
were determined by Mr. Seyfert, the spectral types by Mr. Kuiper. 
The last column contains the linear velocity perpendicularly to the 
line of sight, computed from the angular proper motion and the 
parallax derived from the spectral type under the assumption that 


DM R.A. Dec. mpg m pv Spec. bu Vt 
km/sec 

8 Le eee o537™2 | +71°37’ | 10.30 | 10.34 | Assp | 0736 820 
eee 1:38.9 | +7258 | 10.16 | 10.17 | Assp 132 700 
....... 5:47.1 | +2416 | 10.20 | 10.00 | Fo 830 
15:38.3 | +6708 | 10.46 | 10.46 | F3 .39 520 
42-2667... 15:59.9 | +4232 | 10.18 | 9.94 | F2 630 
18:28.7 | +1303 | 11.00 | 10.63 | A8s 690 
19:08.2 | — 045 9.31 g.o1 | F4 190 
SOOT. 20:05.6 | +42 34 | 10.35 | 9.96 | F3 320 
22:06.7 | +17 36 | 10.00] 9.27 | Fi 0.54 490 


the star lies on the main sequence. If the directions of these motions 
are also considered, it is found that +72°94 would escape from the 
galaxy if it were a main-sequence star.’ The other stars can, on the 
basis of the transverse motions alone, not definitely be said to be 
underluminous. 

On the basis of the spectral features +71°31 and +72°94 are defi- 
nitely intermediate white dwarfs, whereas +13°3683 and +24°1005 
may be intermediate white dwarfs. The remaining stars are prob- 


ably normal stars of high velocity. 
G. P. KurIPer 
C. K. SEYFERT 
YERKES AND MCDONALD OBSERVATORIES 


© Cf. Pub. A.S.P., 47, 280, 1935. The right ascension for the intermediate white 
dwarf +54°2461, given there, should be 20"58™, not 21"58™. 
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SOME NEW SPECTRA OF GALACTIC NEBULAE 


In a recent note’ a preliminary report was made on the spectra 
of some diffuse nebulae. We have now obtained a number of new 
spectrograms (PI. II) which reveal several features of the problem 
of nebular excitation. 

Perhaps the most interesting is the spectrum of the brightest part 
of the “Great Spiral” in Orion, the large diffuse nebula which ex- 
tends through most of the constellation. The slit was centered at 
a 5%45™, 6 + 0°8, and covers 2° in right ascension. The plate shows 
Ha, HB, Hy, Hé, and d 3727 of [Ou]. A defect mars the Ha region 
of the negative, but another plate shows Ha to be very strong and 
confirms all other features.2 The relative weakness of \ 3727 
suggests that the nebula is of low excitation—perhaps even lower 
than the nebulae near y Cygni. The emission character and the 
large extent of this nebula, which is at least 3°5 from any bright star, 
mark it as unusual. If it is assumed that the Orion cluster of B stars 
is the source of luminosity, we can estimate that the linear distance 
of the nebula from the illuminating source is at least 20 parsecs, and 
it may be much greater. 

The next spectrum, that of the nebula near & Persei (NGC 1499), 
is a typical low-excitation emission nebula. The lines visible are 
Ha, HB, Hy, Hb, He, and \ 3727; the nebular lines are weak or ab- 
sent, even though the exciting star is of type O7 and is of intermedi- 
ate luminosity. The spectrum of NGC 1982, the detached knot of 
nebulosity 8’ north of the Trapezium in Orion, is that of a high-ex- 
citation nebula in which Ni and N2 of [O mJ are very strong. The 
visible lines include Ha to Hx, Ni and N2, \ 3727 and A 3869. The 
line \ 3869 of [Ne 111] is strongly concentrated at the center of a con- 
tinuous spectrum, the source of which is doubtful and is in part the 
B8 star near the center of the nebula. A spectrum of the Orion 
nebula at the Trapezium, NGC 1976, does not show the line at 
\ 3869 but is otherwise similar. It is interesting to make a compari- 
son of the North America nebula, excited by a Cygni (Azp), the 
Orion nebula, excited by the Trapezium (O6—-Og) and NGC 1499, 

' Greenstein and Henyey, Ap., J., 86, 620 1937. 


2 We have since obtained a spectrum centered at a 5"50™, 6 —3°. The emission 
features are similar, though the nebula is fainter. 
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excited by € Persei (07). The level of excitation in the North Amer- 
ica nebula is very close to that in the Orion nebula, even though the 
exciting star is of considerably later spectral type; the nebula near 
£ Persei is of much lower excitation even though the illuminating 
star is of early type. This emphasizes the fact that the spectral type 
of the illuminating star is not the sole determining factor in nebular 
excitation. 

The spectrum of NGC 2024, north of and following ¢ Orionis, 
(Bo), shows what is probably a composite spectrum with an un- 
usual change of the relative intensity of the emission lines along the 
slit. Ha is particularly strong in NGC 2024, on the east, but \ 3727 is 
enhanced toward the west, where a faint extension of IC 434 is found. 

We find that the Balmer decrement is subject to surprising varia- 
tions. For example, the nebula IC 2159, surrounding the cluster 
NGC 2175 and excited by an O6 star, shows the hydrogen lines, 
Nt and X\ 3727. A similar spectrum is found for a large nebula near 
S Monocerotis, spectral type O7s. The excitation as indicated by 
the strength of Ni is the same. Nevertheless, Ha is relatively much 
stronger in IC 2159—a fact which is corroborated by the great 
strength of the nebula on direct photographs in the Ha region, 
taken by Struve. 

The spectrum of the Pleiades region was obtained south of and 
preceding Merope. The strong continuous spectrum is partly that of 
Merope and partly of nebular origin. An absorption line spectrum 
extends toward the west for 2°, far beyond the usually photographed 
Pleiades nebulosity. Its character gradually changes with increasing 
distance from Merope, from a spectrum of early type to one of later 
type, since the hydrogen lines Hy and H6 are strongest near Merope, 
while the H and K lines and the G band gain in relative intensity 
at larger distances. The same result is given by measurement of the 
energy distribution in the continuous spectrum which shows a 
change from color class a5 to g5. The early-type spectrum must arise 
from the scattered radiation from Merope, but the late-type spec- 
trum resembles the normal spectrum of the night sky. 

Jesse L. GREENSTEIN‘* 
YERKES OBSERVATORY L. G. HENYEY 


3 Ap. J. 86, 94, 1937. 4 National Research Fellow. 


